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I.   INTRODUCTION 

Recent  developments  in  control  theory  and  practice 
coupled  with  digital  electronics  advances  have  produced 
increased  interest  in  digital  power  plant  modeling.   Although 
conventional  marine  boiler  modeling  is  not  a  new  concept ,  it 
is  still  an  area  of  increased  interest  due  to  the  relative 
lack  of  analytical  data.   This  deficiency  is  due  primarily 
to  the  extreme  environment  that  must  be  endured  by  sampling 
equipment.   Advances  in  data  collecting  and  processing  equip- 
ment have  made  boiler  data  collection  practically  feasible, 
and  in  fact,  nuclear  steam  generator  development  depends 
heavily  on  actual  data.   However,  it  is  not  economically 
feasible  for  a  single  manufacturer  to  retrofit  modern  data 
collection  systems  on  a  standard  marine  boiler.   The  majority 
of  marine  boiler  models  available  are  not  attractive  from 
the  control  engineer's  viewpoint  because  of  one  or  more  of 
the  following: 

a.  The  model  is  over-simplified  to  the  point  of  being  a 
"teapot"  model. 

b.  The  model  is  completely  developed  and  written  in 
Laplace  transforms  and  the  state  space  equations  are 
too  difficult  to  extract. 

c.  The  model  is  over-complicated  to  the  point  of  being 
computationally  inefficient. 
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d .   The  causes  and  modeling  of  "shrink  and  swell" 
phenomena  are  ignored. 

This  model  attempts  to  compromise  between  simplification 
and  complication  while  including  a  theory  and  model  for  shrink 
and  swell.   The  model  follows  the  path  specified  by  Fini  [1] 
restated  below: 

a.  A  general  D-type  marine  boiler  model  is  prepared  in 
Continuous  System  Modeling  Program  (CSMP)  language. 

b.  A  general  FORTRAN  program  is  used  to  prepare  the 
initial  conditions  needed  by  the  CSMP  model. 

c.  The  initial  condition  program  depends  only  on  data 
easily  obtained  from  the  manufacturer's  technical 
manual  and  engineering  handbooks. 


II.   MODEL  CONSIDERATIONS 

A.  BOILER  TYPE 

The  D-type  boiler  was  chosen  for  modeling  because  of 
its  comparative  standardization  among  manufacturers  coupled 
with  its  popularity  in  naval  ship  power  plants.   This 
particular  boiler  is  one  installed  in  the  LHA-1  class  of 
United  States  naval  ships.   Since  D-type  boilers  operate 
under  a  wide  range  of  geometry  and  output  conditions  the 
initial  condition  and  model  programs  are  designed  to  operate 
over  similar  ranges. 

B.  MODEL  NOMENCLATURE 

The  model  nomenclature  differs  from  reference  [1]  because 
of  a  need  to  clarify  the  model  equations  and  afford  readability 
to  the  actual  computer  programs.   The  following  objectives 
were  considered  when  developing  the  nomenclature. 

a.  The  equations  must  be  easily  read,  necessitating 
minimum  variable-definition  referrals  by  the  reader. 

b.  The  notation  used  in  the  model  development  must  be  the 
same  as  that  used  in  the  computer  programs. 

The  resulting  notation  combines  a.  and  b.  and  consists 

of  frequently  used  notations  for  variables  followed  by  a 

subscript  to  indicate  state  point  location  and/or  time  using 

Figure  1. 

Table  1,  a  tabulation  of  general  rules  and  examples, 
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should  suffice  to  allow  the  reader  to  follow  the  model 
development . 

C.  EQUATION  DEVELOPMENT  FORMAT 

The  dynamic  model  and  initial  condition  programs  are 
developed  simultaneously  with  corresponding  equations 
written  in  appropriate  proximity  to  each  other.   A  notation 
of  [IC]  following  an  equation  indicates  that  it  is  an  initial 
condition  equation,  while  a  notation  of  [ME]  indicates  a 
model  equation.   In  some  instances  the  equation  applies  to 
both  programs,  an  indication  of  [IC,  ME]  is  used  to  indicate 
this . 

In  as  many  instances  as  possible  equations  are  developed 
directly  from  elementary  heat  transfer  and  fluid  dynamic 
principles.   This  stems  from  a  desire  to  encourage  more  re- 
search in  this  particular  field  and  to  foster  a  better  under- 
standing of  marine  boiler  dynamics  among  students  of  marine 
engineering.   The  model  was  developed  to  be  read  and  used  by 
students  with  a  minimum  of  additional  research. 

D.  ASSUMPTIONS 

Various  assumptions  are  made  to  facilitate  either  model 
simplification  or  program  tractability .   The  latter  are  not 
always  desirable  and  further  research  or  different  programming 
techniques  may  obviate  the  need  for  them.   The  following 
major  assumptions  are  made: 

1.   The  energy  transferred  to  the  generating  tubes  is 
distributed  uniformly,  i.e.  "even  heating"  exists. 
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2.  There  is  a  thorough  mixing  of  the  fluid  in  the  steam 
drum. 

3.  Steam  generation  occurs  only  in  the  boiling  section 
of  the  screen  and  main  bank  risers. 

4.  All  generating  tubes  and  downcomers  in  a  particular 
circuit  have  the  average  geometry. 

Uniform  energy  transfer  is  assumed  to  facilitate  model 
simplification  and  is  justifiable  since  the  generating 
banks  are  designed  to  optimize  uniform  heat  transfer.   Mixing 
of  the  fluid  in  the  drum  is  a  result  of  the  turbulence  in  the 
drum.   The  third  assumption  results  from  a  lack  of  reliable 
experimental  data  and  the  knowledge  that  D-type  boilers  are 
natural  circulation  boilers  relying  on  a  density  difference 
in  the  circulation  loop  to  promote  flow.   Steam  production 
outside  of  the  generating  banks  would  decrease  the  density 
difference,  thus  inhibiting  natural  circulation.   The  geometry 
assumption  is  necessitated  by  model  simplification. 
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III.   MODEL  DEVELOPMENT 

A .  GENERAL 

The  boiler  is  developed  from  control  volumes  for  which 
the  appropriate  equations  for  energy  transfer  and  storage 
rates  are  written.   The  initial  condition  equations  are 
generated  by  setting  the  rate  (d/dt)  terms  equal  to  zero  and 
solving  for  the  desired  variables. 

The  model  is  developed  with  three  different  elements: 

1.  Furnace-side  heat  transfer 

2.  Water-side  heat  transfer 

3.  Water-side  circulation 

B.  FURNACE-SIDE  HEAT  TRANSFER 

The  gas  flow  path  can  be  followed  on  Figure  1  where  it  is 
represented  by  the  dashed  line.   The  fuel  and  air  mixture 
enters  the  furnace  at  state  point  Q  where  radiation  heat 
transfer  occurs  to  the  screen  riser  tube  metal .   The  flue  gas 
then  leaves  the  furnace  and  exchanges  heat  via  convection  with 
the  superheater  Cstate  point  R-S),  main  bank  riser  tubes 
Cstate  points  S-T),  and  economizer  (state  points  T-U)  in  that 
order. 

At  state  point  Q  the  fuel-to-air  ratio  and  fuel  and  air 
temperatures  are  relatively  constant.   This  allows  the 
sensible  heat  of  the  fuel  and  air  to  be  lumped.   The  mass  flow 
rate  into  the  furnace  is  the  sum  of  the  air  flow  and  fuel  flow 
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The  initial  values  of  fuel  and  air  flow  at  a  specified 
operating  point  are  known.   Therefore,  the  corresponding 
initial  condition  equation  is: 

The  mass  flow  rates  through  the  remaining  sections  of  the 
boiler  are  the  same. 

WQ<3  =  MR*.  -  rviss  -  rn-rr  -  muu.  £me] 

The  total  heat  supplied  to  the  boiler  can  be  written  as 
the  product  of  the  mass  flow  rate  of  fuel  and  the  "lumped" 
fuel  heating  value. 

At  a  specific  operating  point  the  heat  released  to  the  furnace 
is  available  from  the  technical  manual  as  is  the  furnace 
volume.   The  total  heat  input  at  steady  state  is  the  product 
of  these  two. 

QQ0  =  q-ror^    *   fur-vou         trcTj 

This  allows  the  computation  of  the  fuel  heating  value. 
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1 .   Furnace  to  Screen  Risers 

The  heat  transfer  to  the  furnace  screen  tubes  occurs 

primarily  through  radiation,  i.e. 

4-   -r  4 


G)  ■=  <^A(TRR  -Tw) 


which  is  the  Stef an-Boltzmann  law-'-where  o   is  the  Stefan- 
Boltzmann  constant  defined  as: 

Therefore  the  heat  transfer  to  the  furnace  screen  can  be 
written 


where 


Si6m*A  =    <?*  *  k 


At  a  given  operating  point,  technical  manual  values  are 
specified  for  the  furnace  heat  absorption  rate  and  the  area 
of  the  radiant  heat  absorbing  surface  of  the  furnace  screen. 
The  product  of  these  two  yields  the  steady  state  heat  transfer 
rate  which  in  turn  can  be  used  to  compute  SIGMAA. 


This  equation  assumes  that  all  the  radiation  from  the  combus- 
tion flame  strikes  the  screen  tubes,  and  both  the  flame  and 
the  tubes  behave  as  black  surfaces. 
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Sl^MAA  =  Q±^  /((TRR0-h%0.0)*  *  4.0  ~ 
(TW0  4  460.0)  *  *  4.0) 

An  energy  balance  equation  for  the  mass  of  flue  gases  in  the 
furnace  region  can  be  written  in  the  form 

J2.    (  yy\  dp  T)  =  TO*ai  eAQrCN  .Supply d  -  enerovj 
df  Transferred  4o    Screen  4-uJbe.s  — 

or  ^ 

df  *  CjQfL  *  (Tftc-T/Vvt6) 

Note:  TRR-TAMB  is  the  absolute  temperature.   At  steady  state 
the  rate  terms  are  equal  to  zero.   Values  of  total  heat 
supplied  to  the  boiler,  heat  transfer  rate  to  the  screen 
tubes,  total  flue  gas  flow  rate,  and  flue  gas  temperature 
are  available  from  previous  equations.   This  allows  the  specific 
heat  of  the  flue  gas  in  the  furnace  region  to  be  computed. 

&3  fc  -  (QG^  -  qi0  )  /  (ttfc/<$  *  (TQfiUfi  -TAfAfc))         [it] 

During   a   transient,    the   density   of   the    flue   gas    in   the   furnace 
and   the   specific    heat    of   the    furnace    flue   gas   can   be   regarded 
as   constants   near   a   steady   state   operating  point.    The   dynamic 
equation   can   be   written   as 
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where 
Thus, 


or  in  CSMP-III  language 

The  value  for  TRRO  is  available  from  the  technical  manual 
2 .   Flue  Gas  to  Superheater 

For  the  remaining  tube  banks  heat  transfer  occurs 
primarily  through  convection.  The  energy  given  up  by  the 
flue  gas  in  flowing  over  the  tube  surfaces  is  represented 
by  the  general  equation 


energy   =  mass  flow   *   specific 
given  up   rate  of       heat  of 
flue  gas      flue  gas 

*  change  in  temperature  of  the  gas 


1.  e 


r;  rv\  Cp  M~" 
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At  steady  state,  the  energy  transferred  to  the  superheater 
tubes  from  the  gas  must  equal  the  energy  transferred  from 
the  tubes  to  the  steam. 

The  steady  state  values  for  superheater  inlet  and  outlet 
steam  temperature  and  pressure  are  specified  by  the  technical 
manual  along  with  the  superheated  steam  mass  flow  rate.   This 
allows  the  direct  computation  of  Q40  and  Q30. 

G^0  =  MV\n0   %  (muiv)^  -  HMM0)      [IC] 

HMMO  is  the  enthalpy  of  saturated  steam  corresponding  to  drum 
pressure,  computed  with  curve  fitted  equations  [2].   HNNO  is 
available  directly  from  the  technical  manual.   The  specific 
heat  of  the  flue  gas  in  the  superheater  region  can  now  be 
computed. 
COS   -  Q30  /  (P\R«$   *  (TRRJ0  -TSS0Y)   tCGji 

Thus  the  dynamic  equation  for  heat  transferred  from  the  flue 
gas  to  the  tube  metal  is  written: 

The  heat  transferred  to  the  tubes  via  convection  may  also  be 
represented  by  the  Grimson  correlation  [3]  which  states  that 
for  cross  flow  over  tubes; 
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^  K     %     J 

where 

U.00  =  Veloc.'ty 


Since  h  can  also  be  written 


r* 


<\ 


then  / 


also 


leading  to  the  general  equation 


K, 


^^hrfx-)  ^ 
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Around  a  particular  steady  state  operating  point  the  physical 
properties  can  be  regarded  as  a  constant  and  the  physical 
dimensions  are  constant,  i.e. 

AC   (—A  V  ft-'^     =  Coosiarff 

which   implies 


q    ^    fcons+an+)  #•  A  "  /ST 


For  flow  over  tube  banks  with  differing  rows  of  tubes  and 
differing  tube  geometries  the  constant  (n)  varies  little  and 
the  average  is  approximately  .6  [3j.   Therefore,  in  the  case 
of  the  superheater  and  main  bank  riser  tubes,  the  correlation 
is  used  in  the  form 

or  for  the  superheater 

Q3  -  KRS  %  MRfiL  ^  *  .6  *.  (TR.£  -TWW) 

where  TRS  is  the  average  flue  gas  temperature 

Since  the  steady  state  values  for  Q30,  MRRO ,  TRRO ,  TSSO,  and 
TWO  are  given  by  the  technical  manual  the  constant  KRS  can 
be  calculated. 
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To  provide  the  dynamic  solution  of  superheater  flue 
gas  outlet  temperature,  the  two  transient  equations  for  the 
flue  gas  to  superheater  tube  metal  heat  transfer  are  equated 
resulting  in 

TR.R.  -  TSS 


where 


trs  -  (jfcje.  -v-tss^/^.0 


Solving  the  above  for  TSS  gives 


tss  =.  (Vfcfc.  ^  (Pttit-l)+z  ^twu;)/Cphxi+1)    \W\e] 


where 


Phti  -  z,  *  fv\*tfL*  *  .M  *  cj^s/k^s        Hm£1 

3 .   Flue  Gas  to  Main  Bank  Risers 

The  equations  for  the  flue  gas  to  tube  metal  energy 
transfer  in  the  main  bank  are  derived  similarly.   For  the 
equation  involving  flue  gas  specific  heat; 

The  Grimson  correlation  equation  for  the  main  bank  risers  is 

qr-a  ksv  £  mss  *.  Xs  fl.(b  ik  (tst  -tyv) 


where  . 

Tsr  =  (tss  +ttt)  /?-9$ 

By  equating  these  two  equations  as  was  done  with  the  super- 
heater the  flue  gas  temperature  leaving  the  main  bank  can 
be  written 

ttt-(tss  *  (pHrz-l)-h^^Tvv)/(PHr^^-l)  [we] 
phiz  -  &*  hss  *  *  0.4  ^  cst/kst     [meJ 


where 


The  steady  state  heat  transfer  is  computed  differently  than 
that  for  the  superheater.   Since  the  steady  state  heat  trans- 
fer to  the  screen  bank  has  been  computed,  the  energy 
transferred  to  the  main  bank  is  the  difference  between  the 
total  energy  transferred  and  the   desuperheater  and  screen 
bank  energy  transfer.   At  steady  state  the  mass  flow  rate  of 
steam  out  of  the  boiler  is  equal  to  the  mass  flow  rate  of 
feed  in  the  boiler  and  the  total  energy  transferred  is  the 
mass  flow  rate  multiplied  by  the  enthalpy  change  between  the 
outgoing  steam  and  incoming  feedwater.   Thus,  the  steady 
state  energy  transfer  equation  for  the  main  bank  is 

Q990  is  computed  in  the  water  side  heat  transfer  section. 
Since  Q50  =  Q60  at  steady  state  the  heat  transfer 
coefficient  KST  and  the  specific  heat  CST  can  be  computed. 

Ksr  «  qs-0  /  (mssjzS  *  *  0.<*  *  (j^jp  -rr/0))  trc] 
Csr  =g!s#/(mss0*  (rsstf-rrrfi))     [re] 
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4 .   Flue  Gas  to  Economizer 

The  economizer  tubes  differ  from  the  generating  and 
superheater  tubes  in  that  they  are  finned.   A  reasonable 
finned  tube  heat  transfer  correlation  was  developed  by 
Weirmann,  et  al .  [4] 

'111 

Co  Sr*a* 


J    -   -^ 


G     =  maximum  mass  flow  rate 
max 


j  e  correlation  constant 
C   e  specific  heat  of  flue  gas 
Since      n  ~zl      Q  the  above 

correlation  can  be  written 

Q  =  j  A  A.T  tf>  6.  wax 

pr*/3 


G    is  a  function  of  mass  flow  rate  into  the  tube  bank  and 
max 

tube  bank  geometry.   Similar  to  the  case  of  the  Grimson  coeffi- 
cient used  for  the  superheater  and  main  bank  tubes  the  physical 
and  geometric  properties  are  considered  constant  leading  to 
the  use  of  the  correlation  in  the  form 


^ 


--    CmAT 
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Therefore,  for  the  economizer  the  equation  is 

Q7    ~  KTu.  *  mxt  *  (ttu.  -txx) 

where 

TTUL    =   TTT  -1-TU.U. 

Paralleling  the  superheater  and  main  bank,  the  heat  transfer 
may  also  be  written 

Q7=  mtT  *  CTu  ^  (TTT-  T"ulul)  [tv\eQ 

The  steady  state  equation  development  follows  that  of  the 
superheater  and  main  bank  sections  and  is  not  repeated  here. 
The  equations  are: 

Gaj>  -  Q*<?  LreJ 

cm  -  Qn4>l(mr<t>  #.(tttj6  -  tllu.^))  [k3 
Ktu.  =  Q.i<pl(nvr0  *(ttu.0  -  T*a0))   r_rcQ 

TTU0  -  (T-TT-0  -*-  TU.U0 )  /z.  0  jjnfj 

Knowing  the  values  for  KTU  and  CTU  developed  by  the  initial 
condition  program  allows  the  solution  of  the  flue  gas  outlet 
temperature  in  the  transient  analysis.   Thus, 

Imej     Tau  =•  (rrr  5K.  (pHr3  -Lp)+Z.0*k  tm)/ 

(PKT3  +1.0) 
[me]       PH-I3   =  2.0*.    (LTU./KTLL 
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C.   WATER  SIDE  HEAT  TRANSFER 
1 .   Risers 

The  heat  supplied  to  the  riser  banks  is  used  both  in 
the  boiling   and  nonboiling  length  of  the  tube;  however  the 
assumption  is  made  here  that  the  heat  transfer  is  all  under 
fully  developed  flow  boiling  conditions.   Since  the  nonboiling 
length  of  the  risers  is  relatively  small,  little  error  is 
produced.   For  fully  developed  nucleate  flow  boiling,  Levy 
[5]  suggests  the  equation 

\.T$X  *  )0fe  sec- tT 

By  lumping  the  constant  area  and  denominator  the  riser  equations 
can  be  written  in  the  forms 

Therefore,  the  two  riser  equations  are 

Q2.  -  kV  *  PSAT*  *  (4.<?S/M)  *. (yW-Tshr)*.  *•  3.0   LmeJ 

and 

Qfc  ^  KY*.  few  **(n.0/3.0) *(r/Y  -Tw)*  *3.0  M 

for  the  screen  and  main  bank  respectively.   A  simple  energy 
balance  on  the  tube  metal  yields  the  metal  temperatures,  i.e., 
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The  equations  for  screen  and  main  bank  metal  temperatures  are: 
TVV/  -  II^T^L^Tl/l/^OT/l/)  ^£] 

orvv  -  (qs*-  qc*)  /  Classy  *  cpy)  ]W] 


The  heat  transfer  coefficient  KY  can  be  computed 
directly  from  tube  geometrical  data  available  in  the 
technical  manual; 

Since  Q10  =  Q20  at  steady  state,  the  dynamic  equation  for 
Q20  can  be  solved  in  a  steady  state  situation  for  KV, 


[ 


id    kv«  qi^/(Ps^T**  (4/3)*  (Vw0-Tspfr)#*3.0) 


Q60  was  computed  in  Section  III-B.3.  This  allows 
the  calculation  of  the  initial  main  bank  riser  tube  metal 
temperature, 


ic]  tw0 = fa<o<f>/(xiK  es*r**  6/3)))  *  *.  (1/3)  +■  t^t 


In  the  case  of  the  screen  risers,  the  initial  tube 
metal  temperature  is  available  from  the  technical  manual  and 
the  initial  heat  transfer  rate  has  been  previously  calculated 
This  allows  the  calculation  of  the  screen  riser  steamside 
heat  transfer  coefficient, 


icj     K^  *  Q10/(RSAT  ¥  *  (M .0/3.0)*  (tW0- TSffi)**  3.0) 
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2 .   Desuperheater 

The  Dittus-Boelter  correlation  [6]  is  used  to  compute 
the  heat  transfer  rate  from  the  steam  to  the  desuperheater 
tubes.   This  correlation  in  its  basic  format  is 

OX. 


NJu.  -  .^Z3  R^'Prr 


where 


Nu  =>  Nusselt  number 

Re  =>  Reynolds  number 

Pr  =>  Prandtl  number 

n  =>  .3  (for  cooling  fluid) 

When  the  appropriate  variables  and  constants  are  substituted 
for  the  dimensionless  numbers,  the  Dittus-Boelter  correlation 
can  be  written 

where  the  constant  K  is  defined  as; 

K  -  -0ZZ  Ik   Uermal  concLAitfty   %  taai"  W<sTer  area 


VlT; 


and  AT  is  the  log-mean-temperature  difference  (LMTD) 
Thus,  for  the  desuperheater,  the  equation  is 
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q°\    -     KlOP    *.   f^KJNJ  &  *   .S  *.    LMTONiP 

KNP    can   be   computed   directly.      Thus, 

KMP  ^  (( . <#ZZ  &  THtojOfj) /oMp) *i(i\.<fi/(PZ*  DlOP 

#;    AR.EAOS  ^Q 

The  energy  given  up  by  the  steam  to  the  tubes  may 
also  be  written  in  terms  of  steam  specific  heat,  flow  rate, 
and  temperature  difference. 

£^  ~  £^P  *  MKJrO^(xfOiO  -Tpp)  I^eH 

Similarly  since  the  heat  transfer  to  the  water  in 
the  drum  is  totally  by  convection; 

For  initial  condition  calculations  THE  is  considered 
to  be  equal  to  TSAT.  Therefore 

The  desuperheater  outlet  temperature  is  solved  by 
equating  the  two  dynamic  equations  for  Q9  and  solving  for 
the  desired  outlet  temperature. 

TPP  =  (TMNJ  -  Ttt)/ (e*P (MP /  (<LW(°  *  MWU**  fi.Z))) 
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The  instantaneous  desuperheater  tube  metal  tempera- 
ture is  obtained  by  an  energy  balance  on  the  tube  metal 
similar  to  that  for  the  riser  tube  metal. 


r^ass  *i  dp  K  d~T  -  q.  -  q  • 


Solving  the  appropriate  equation  for  dT/dt  results  in 

or**  =  (§°i-G{<W)/(m$s*  x  ops)     j>e] 

The  specific  heat  of  the  steam  in  the  desuperheater 
(CNP)  may  be  calculated  using  the  known  steady  state  values 
for  desuperheated  steam  flow  rate  and  the  desuperheater  inlet 
and  outlet  temperatures  and  pressures. 

ttjp  •=  qS^/(m(Jio{zS^.(tmm^-tpp0))  [re} 

With  the  value  of  Q90  computed  the  initial  log-mean- 
temperature  difference  can  be  calculated. 

The  log-mean-temperature  difference  is  a  function  of 
steam  temperature  in,  tube  metal  temperature,  and  steam 
temperature  out;  therefore  the  initial  tube  metal  temperature 
may  be  calculated. 
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3 .   Economizer 

Paralleling  the  desuperheater  development  the  Dittus- 
Boelter  correlation  is  used  to  relate  the  heat  transfer  from 
the  tube  metal  to  the  feedwater.   Thus, 

the  only  difference  being  that  n  is  now  .4  vice  .3  because  the 
fluid  is  now  being  heated.   The  appropriate  constants  are 
again  lumped  yielding 

Q?  -  KX  &    Wl\t\  *  ^  0.?  3fc  UVAT0A6 

An  energy  balance  on  the  water  gives  the  heat  absorbed  in  two 
other  forms  ; 

q?  ^  roAA  *.(h<56  -haa) 

q%  -  maa  *.  ca€  *.  (Tise  -taa)  0*e1 


At  the  specified  operating  point  the  inlet  and  outlet  condi 
tions  of  the  economizer  as  well  as  the  mass  flow  rate  are 
given,  allowing  the  computation  of  the  steady  state  heat 
transfer. 

QZ0  ■=   MAA0  *  (t+S6^  -HAA0)    I^c] 
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The  heat  transfer  coefficient  KX  is  computed  directly. 

kx=  ((.^z3*thcon;a)/oa6)#  O4.0/(pr*oA8a< 

^  AfceAEC  YrCl 

This  permits  the  solution  of  the  steady  state  log-mean-tem- 
perature difference. 

uwro  AS  a  q  §-0  / (XX  *  mAA^  *  *  $  •  ?)     &<£! 

Following  the  identical  path  delineated  for  the  desuperheater , 
the  tube  metal  initial  temperature  can  be  computed. 

E>(POBL  -  E)CP((TT8(S^-TAAt^  )  /  L^TO (\&  l^c] 


An  energy  balance  on  the  economizer  tube  metal  yields 
its  instantaneous  temperature. 

dtxx  -  (Q?-q7)/(^Ass)C  *  cpa)      l>e] 

T"KK    =  ItJTe.KU^T/X^DTxx)  £.™ej 

By  equating  the  formulations  for  Q8  involving  LMTD 
and  specific  heat,  the  economizer  outlet  temperature  may'  be 
calculated . 

T66^  (TAA  -TX*)/(e*P  (ka/(cA6*MAA**  <j>-z))) 
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The  specific  heat  of  the  feedwater  in  the  economizer 
is  calculated  at  steady  state  using  the  computed  value  of 
steady  state  heat  transfer  and  the  given  mass  flow  rates  and 
inlet  and  outlet  temperatures. 

4 .   Superheater 

The  development  of  the  superheater  equations  follows 
that  of  the  desuperheater  and  economizer  and  for  that  reason 
the  development  will  not  be  repeated.   The  equations  are 
listed  below. 


The  initial  superheater  log-mean-temperature  difference 
is  computed  directly  allowing  the  subsequent  calculation  of 
the  tube-metal-to-steam  heat  transfer  coefficient  used  in  the 
Dittus-Boelter  equation. 

Lw\TOrv\K)  -  (TUM0-TMlv\0V(ALoa  ((TWW0-Twrt0)  / 

(tvv/w^-tmn;^)))  \xc] 
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An  energy  balance  on  the  tube  metal  yields  the  instantaneous 
tube  metal  temperature,  the  initial  tube  metal  temperature 
being  available  from  the  technical  manual. 

DrwW  -  (q3_q4)  /(^assn*/^  cpW)       \>e] 

TVs/Vs/       -     lUTSftU    (tW\a/0^  OTVv/U/)  Xj^^l 

The   superheater   outlet    temperature   can   now  be   calculated. 

tnjnJ  -.  (rrA^-TWiA;)/(EXP  (kW  /(cmJ  *  niUM*  * 

D.   WATER-SIDE  CIRCULATION 
1 .   General 

The  water/steam  side  circulation  equations  are  by  far 
the  most  difficult  to  visualize  and  codify.   These  equations 
must  be  accurate  in  order  to  predict  phenomena  such  as  shrink 
and  swell  while  appropriate  assumptions  and  simplifications 
must  be  made  in  order  to  make  the  equations  tractable. 

As  can  be  seen  in  Figure  1 ,  the  f eedwater  enters  the 
economizer  at  state  point  A,  passes  through  the  economizer 
to  the  steam  drum  where  it  becomes  part  of  the  water  volume. 
The  liquid  leaves  the  drum  via  the  downcomers,  both  main  bank 
and  screen  at  state  points  G  and  C  respectively.   The  main 
bank  downcomer  delivers  its  liquid  to  the  main  bank  risers 
via  the  water  drum,  and  the  main  bank  risers  deliver  the  fluid 
to  the  steam  separators.   The  flow  through  the  screen  risers 
is  the  same  with  the  exception  that  there  is  no  water  drum  in 
this  circuit.   The  steam  separators  separate  the  majority  of 
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the  water  from  the  steam  leaving  the  drum  at  state  point  M 
and  the  majority  of  the  steam  from  the  water  being  discharged 
back  into  the  drum  liquid.   The  steam,  with  a  negligible 
amount  of  water  passes  through  the  superheater  via  state 
point  M-N.   At  the  outlet  of  the  superheater  the  steam  leaving 
at  path  III  is  used  for  "main"  steam;  that  leaving  via  path 
II  travels  to  the  desuperheater  where  it  is  cooled,  then 
leaves  via  state  point  P  and  is  used  for  auxiliary  steam. 

The  water  leaving  the  separators  with  a  small  amount 
of  entrained  vapor  is  mixed  with  the  incoming  feed  water  and 
forms  a  "foaming"  vapor/liquid  mass  in  the  bottom  half  of 
the  drum.   The  liquid  from  this  mass  leaves  via  the  downcomers 
and  is  circulated  through  the  loop. 
2 .   Downcomer  Pressure  Drop 

A  closer  look  at  the  downcomer  flow  is  necessary  to 
justify  assumptions  that  are  made  in  the  development  of  the 
pressure-drop  equations. 

Circulation  ratio  is  defined  as  the  ratio  of  the  total 
weight  of  steam  liberated  to  the  drum  [7].   For  a  600  pound 
marine  boiler  this  circulation  ratio  is  on  the  order  of  20:1 
{7  and  8] .   This  implies  that  for  every  21  pounds  of  liquid 
flowing  down  the  downcomers,  20  pounds  of  it  has  already 
traveled  up  the  riser  and  is  at  saturation  temperature.   There- 
fore, assuming  the  downcomers  are  perfectly  insulated,  it  is 
reasonable  to  assume  that  the  downcomer  liquid  is  at  or  very 
near  saturation  temperature. 
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The  momentum  equation  for  the  downcomers  may  be 
written : 

pressure  at  top  -  pressure  at  bottom 
=   frictional  loss  -  gravitational  head 
+  entrance  loss    +  bend  loss 
+  exit  loss       +  inertia  force 
The  inertia  force  term  may  be  considered  negligible  [9] . 
This  is  the  quasistatic  approximation  which  basically  states 
that  pressure  waves  move  much  more  rapidly  through  the  system 
than  the  important  time  constants.   This  is  mathematically 
equivalent  to  the  elimination  of  a  large  negative  eigenvalue. 
This  quasistatic  approximation  is  only  good  for  on-line 
transients  and  does  not  apply  for  extremely  large  discon- 
tinuities which  would  result  in  the  boiler  being  taken  out  of 
operation,  (e.g.,  a  ruptured  tube).   Therefore  the  downcomer 
momentum  equation  can  be  written: 

pressure  at  top  -  pressure  at  bottom 
=  (friction  factor  *  downcomer  length 
t  downcomer  diameter  +  entrance  factor 

+  bend  factor      +  exit  factor) 

2 
*  (downcomer  mass  flow  rate)  i    (2  *  (downcomer 

2 
cross  sectional  area)   *  density  of  fluid 

in  downcomer  *  g  )  -  density  of  fluid  in 

downcomer  *  gravitational  acceleration  * 

height  of  downcomer  r    g 
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In    the   model    notation 


PtPD  -  (PCO^  Ldo/oeD    -)-    EtOTf^O)  -V- SEnJOCO 

P6PK    *(f=6H   *  LGtt/06tt  -HEn^GB  +  GEkIdSB 

+  EXlT6rt)*  l*AS$  **  2.0/(2.0  sK  ASH *#Z.0 


for  the  screen  and  main  bank  downcomers  respectively.   FCD  and 
FGH  are  friction  factors  of  the  form 

4  *   4-  /(l.74  -fc  LOS  («/KS» 

where  R  is  the  pipe  radius  and  KS  is  the  relative  sand  rough- 
ness 111]  . 

fcd  -  1.0  /(.l.rW  -z- ?5  #  alos i<5z5  (kscd))    IjCjMe] 

F6H  =  1.0/(l.-H  -2.JZ5  *  ALO6I0(ks<5H))  &JL,MeJ 

3 .   Riser  Pressure  Drop 

The  momentum  equation  for  the  riser  boiling  section 
must  take  two  phase  flow  effects  into  account  because  the 
flow  in  the  boiling  section  is  not  homogeneous.  Thus,  there  is 
a  relative  velocity  between  the  liquid  and  vapor  phases  here. 
The  steam  separators  are  included  in  the  riser  length.  However, 
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the  effective  length  of  the  separators  and  thus  the  pressure 
drop  are  considered  negligible  because  of  a  general  lack  of 
information  concerning  these  items. 

The  pressure  drop  due  to  single  phase  or  homogeneous 
flow  in  the  nonboiling  riser  section  can  be  written 

APOT^„  =  AP     ,    .  .    +  AP-  .  .  .    +  AP 

SPF     acceleration     friction     gravity 

Experiments  conducted  by  Babcock  and  Wilcox  [10] 
indicate  that  these  homogeneous  flow  pressure  drops  may  be 
modified  to  give  the  appropriate  two  phase  flow  pressure  drops 
using  correction  factors  that  are  functions  of  outlet  quality 
and  operating  pressure.   The  two  phase  flow  pressure  drop  can 
be  written; 


AP^^  =  AP     ,    .  .    *  r     , 

TPF     acceleration    acceleration 

+  AP  *  r 

friction        friction 

+  AP    .        *  r 

gravity         gravity 


where  the  r  terms  are  two  phase  flow  correction  factors 
available  by  using  curve  fits  of  the  data  from  reference  [10] 
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The  form  of  these  r  terms  become 


4- .  £734 


c*a 


ttKlCK.  -  -34.  ^?£&  *  XLL0  *  *  *.j>  +25.  7J*W  *  XLL^S     IjDtfl 


Therefore,  in  the  boiling  region  of  the  riser 
The  total  pressure  drop  across  the  length  of  the  riser  is; 


AP  =  APSpF  +  APTpF 


In  model  notation  the  equations  are 
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fcuooo  ^  Aoe  #  *  2.^)  +-O-I.05  3<  fde.  *  l£& 

#:   MOD  #  *  2.  j2$  )  /(<££  *.  ODE.  *.  (RHODD  +■ 

RHoee)  >fc  Aoe  #  *.  2.  fi)  \  (s  x^oe.  *.  i^^o 

4- KHOEe))/ (66*2.0)  -i-  (^00**2-0  * 

fcACLE.)/C<£C  *  I^HOEE.^  AEF  *  *  2.0)    +- 
fe.0   *   FEF  *  LEF^  mDD  *  *  Z^,*.  tfFKJC£.)/ 

C$C   *   OEF  *  fcHcSE.  Hi  AEF  #  ^  2.^)  +- 

(&  *  £EF  ^  J^HOEE  yk  /K€A\/e)  /  £C 


FKL  *    LKU*  W;T#*3.0  #   flFRlCfc.)/ 
fee  *  DKL  *   tfMoKK  *  /VRU  *  -*  *.0)  +- 
(6  *  £ku>£  £I4oKK3<  ft$£M\C)/6<L 
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The  friction  factors  are  in  the  same  form  as  those 
for  the  downcomers 
FOE-  l.jty(i.7N  -*.^  *  ALOSIJZS  (KSOE))      M 

Fox=  1.0 /(l/H  -^.^  *  aud£i$zS(ks:jiO>)       \rcj 

At  a  specified  steady  state  operating  point  the 
pressure  drop  across  the  downcomers  must  equal  the  pressure 
drop  across  the  risers  and  the  downcomer  flow  rate  must  equal 
the  riser  flow  rate.  Therefore  by  equating  the  appropriate 
pressure  drop  equations  and  solving  the  resultant  relation  for  the 
flow  rates,  the  initial  flow  rates  may  be  computed  as 

+-fttoEE^)/2.0)-G  *  £BF0**HOEE0*fc 

GE.0DCO  -V-  HX  ITCO)/  (2^  *  ACO  *  *  Z-0 

*  *HoC/>0)4-((ft>to©E0-RttO0O0V 
(RHOE60  X  £HCOO0  *  AOF  *  *.  Z.0  V)  +• 
(4.^3*   FDE    *  LDE£  *  2.0)/Cz.fl  XODF 

jfc  (khoee^  +  khodo^)  tK  adf  #  #  2,0)  -H 

KACL&/(KttoEE0  #  ADF*  *  2.0)  +  6l-0 

*  FEF  *  IEF0  *  ^F^|CE.)/(z.jr5  *  DD^ 

*  R^E£jZ$  *  ADF  *  #  *.0ft)  #  *  0.5-      [rc] 
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WU.00  =   (XftH06VV^*S*ft6  4-<S*2TK^*((KtfOJ^ 

ft£GAA/K0/((F6H  *  L<SW/06itf +■  Ekin^H  -V- 
6£nJ06H  -t-Exrr$H)/6u0-5K  A^H*:*2.0 

(|RHCKK{Z5  *  KHOJT0  *  atl  *  *  2.0})  +- 


4 .   Riser  Continuity 

The  relationship  between  the  riser  inlet  and  outlet 
mass  flow  rates  is  written  in  terms  of  the  continuity 
equation  for  one  dimensional  unsteady  flow. 


In  model  notation,  this  becomes 

tALL  -  m^X-  DRBOTL  *   VOLTU  If^] 

»*\FF    =.  rt£D~  DRH-OOP  X  VOUDF  tj^Cj 

The  numerical  differentiation  technique  used  by  CSMP- 

III  is  highly  suspect,  as  are  other  techniques.  An  "averaging" 

system  is  used  in  the  actual  dynamic  model.   In  addition,  the 

flow  rate  down  the  downcomer  and  into  the  riser  is  held  con- 
stant for  the  open  loop  boiler  model. 
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5 .   Riser  Quality  and  Density 

The  average  density  in  the  risers  must  be  solved  for 
explicitly.   Linearly  varying  quality  along  the  riser  tube 
length  follows  directly  from  the  assumption  of  uniform  heat 
flux  along  the  riser  tube  length.   The  average  density  along 
the  tube  length  is  the  sum  of  the  total  change  in  density  and 
the  density  entering  divided  by  the  riser  tube  length  II]. 
Assuming  the  density  varies  only  in  the  boiling  length  of  the 
riser,  the  total  change  in  density  can  be  written; 

The  average  density  is 


where     p_  =  density  of  saturated  water 
L  =  total  tube  length 


Since  p(il)  varies  linearly  from  entering  to  exiting,  p(l)    can 
be  written  in  the  form 

L-s  ~ 


where    y   =  specific  volume  enterinj 

X    =  quality  at  riser  outlet 
out    M      ■ 
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LR  =  boiling  length 

LND  =  nonboiling  length 

V    =  change  in  specific  volume  from  saturated 
g    liquid  to  saturated  vapor 

The  integral  of  the  explicit  equation  for  p(£)  may 

be  solved  after  rearrangement  to 


ott)    - 


JL 


1-6      '     Le 


which  is  of  the  form 


\ 


s 


av 


This  yields  the  average  density, 


+ 


ft 


L 


10(3 


In  model  notation  the  formula  is  written 

RHODF0  *(±.0/uOF)*  (lef0/((xff0)#  VfiO 

*.  ALO$(((XFF0)/\fl=)*>/F§   -hl.0^ 


43 


To  permit  calculation  of  the  quality  term  the  energy 
balance  on  the  liquid  in  the  riser  tubes  is  evaluated.   Thus, 


rate  of  change  =  flow  energy  -  flow  energy 
of  riser  energy      in  out 


+  thermal  energy  in 


ax 

If  the  enthalpy  term  is  treated  as  the  average  enthalpy  in 
the  riser  the  equation  may  be  written 

In  model  notation  this  is  written  in  the  form: 

Ohxll  -(was  ^(htt-hf-xll.x^fgx/z.^)+-  Q<q  -^ll 

-rv\FF  *  XFF  #  UFS/z.^)/(fcrtOOF#V0LEF)   [>EJ 
H*FF  =  liOT^euCv^^FF^^HXFF)  Jjvve^ 

XFF  -  (u-XFF-  Hf)^^.^/hF6  ^v\e] 

The  initial  conditions,  HXLLO  and  HXFFO  are  calculated  in  the 
initial  condition  program. 

The  steady  state  equations  for  quality  are  derived 
differently.  The  energy  transfer  rate  to  the  riser  liquid 
may  be  written  in  the  form 
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where 


»  ^£  +  ^o J-  ka   -  k;o 


Solving  for  X    yields 

out  J 


*o^-  6$+*(kf-^))/fa  k%) 


or 


XFF0  -  (ql^  f  f*FF0  *  (HF  -H0O^))/(mff^  *  HF6v)      [xc] 
A  LL^>  =  (^5-fS  +■  rv\LL(^  *  (>F  -  HTT<p ))  /  (*LL#  *  HF<k)    JnQ 

6 .   Riser  Boiling  Boundary  Location 

The  nonboiling  length  of  riser  tube  is  the  product 
of  the  riser  tube  length  and  sensible  to  total  heat  ratio  [1] 

At 

The  sensible  to  total  heat  ratio  may  be  expressed  as  the 
enthalpy  change  in  raising  the  water  to  saturated  conditions 
divided  by  the  total  enthalpy  change,  i.e. 

The  equations  for  nonboiling  length  are 

■V-  XFF jZ$  #.  RF$)  -  HDD0  )  ia] 

+-XU-0  ^    HF$)  ~  HvT3#)  \T(fl 
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The  Fortran  function  AMAX1  is  used  here  because  the 
above  equations  are  used  in  an  iterative  loop  in  the  initial 
condition  program  and  there  is  a  possibility  during  iterations 
of  reaching  a  situation  where  the  enthalpy  entering  is  greater 
than  saturation  enthalpy.   The  final  initial  condition  solu- 
tion prevents  this.   The  AMAX1  function  is  not  used  for  the 
dynamic  model.   Thus, 
toe  -    LDF  *  (hf  -Hob)  /((  HF  -hXFF  *  HFS)-  HDo)  \W»b] 

UK  =    LTL  *  (HF  -  HJT)/(C  HF  -VALL  *  HF($)  -HTT)    [we] 

The  boiling  volume  is  derived  directly  from  the  solution  for 
the  nonboiling  length. 

,   . ,  .      ,       .  .  ,    ,     ,,.  boiling  length    „ 

boiling  volume  =  total  volume  *  5 p —    Hence, 

total  length 

VCLEF       =  VOLDF     ^     UEF/LDF  |>e] 

VOlXL.        -  VOLT\.     #       LKl/lTU  JViE] 

7 .   Steam  Drum  Liquid  Mass  and  Energy  Balance 

The  rate  of  change  of  liquid  mass  within  the  drum  is 
equal  to  the  sum  of  the  mass  flow  rates  entering  and  leaving 
the  drum.   The  liquid  mass  flowing  into  the  drum  is  the 
saturated  liquid  from  the  risers,  the  liquid  from  steam  con- 
densation in  the  drum,  and  the  incoming  feedwater.   The  liquid 
leaving  is  that  leaving  via  the  downcomers.   The  dynamic  model 
equation  is 
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The  instantaneous  mass  of  liquid  in  the  drum  is 

where 

The  drum  energy  balance  is  derived  similarly  (  i.e., 

DDwOHL  =wtu  #  (i^J-Xul)^  VVF  +r^\FP^  (i.0-XF=f) 
#  bVP  +  nCo*JO  *  HF6  -V-rvvS6  #.  H66  -  ^^C 


The  instantaneous  drum  energy  is 

The  initial  drum  energy  is  the  product  of  the  initial 
drum  liquid  mass  and  initial  drum  enthalpy. 

The  initial  drum  liquid  enthalpy  is 

OH    -  OmDHL  /omfvs z. 

In  the  dynamic  model,  drum  liquid  enthalpy  is  considered  the 
enthalpy  of  saturated  liquid. 
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The  steam  condensation  rate  equation  is  based  on  the 
difference  between  the  pressure  and  temperature  of  the  steam 
and  that  of  the  liquid  [12] 

-^S/VT/(TSAT4-H(b^.^)-3K^^.^)i-  $ZS<ot 


The  water  level  in  a  marine  boiler  is  generally  con- 
sidered in  reference  to  a  mid-drum  zero,  that  is,  a  water 
level  of  plus  one  inch  implies  the  water  level  is  one  inch 
above  the  drum  centerline.   For  computational  efficiency  the 
assumption  is  made  that  for  small  changes  in  water  level 
around  the  midpoint  of  the  drum  the  water  surface  area  remains 
constant.   This  allows  a  simplified  level  equation,  i.e., 

DMOV  is  the  equation  for  the  total  volume  of  "liquid"  in  the 
drum.   Recalling  that  a  small  percentage  of  steam  leaving  the 
separators  is  entrained  in  the  liquid,  the  rate  equation  for 
"liquid"  volume  in  the  drum  is 

DDrOOV  ■=  ((wFF  +  WLL-WGe^  VF  +  PCu.  IK  (mll-H*)Ff) 
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where  DMDVO  is  the  volume  occupied  at  time  zero  which  is  half 

the  drum  volume. 

8 .   Circulation  System  Initial  Condition  Iteration 
Procedure 

The  initial  conditions  for  the  circulation  system 
are  found  by  flow  rate  balancing  in  the  downcomer-riser  flow 
loops.   An  initial  guess  of  flow  rate  is  made  using  the  assumed 
riser  exit  quality  of  .05  percent  which  is  reasonable  for  this 
type  of  boiler  {7].   Coupling  this  assumption  with  the  assump- 
tion that  the  percent  carryunder  is  zero,  an  initial  approxi- 
mation of  the  flow  rate  can  be  determined  using  the  figures 
for  initial  heat  transfer  rate  to  the  risers  calculated  in 
Section  III-A.   Hence, 

heat  transfer  rate 


mass  flow  rate  = 


assumed  quality  *  latent  heat  of 

vaporization 


A  first  approximation  of  downcomer  enthalpy  for  both 

banks  can  be  calculated.   This  downcomer  enthalpy  is  assumed 

to  be  the  same  for  all  downcomers .   Using  an  energy  balance 

on  the  liquid  in  the  steam  drum 
MFF^  s  QA^/(XASur*lE.*HF6i)         &CI 

The  screen  riser  inlet  enthalpy  is  assumed  to  be  equal  to 
that  of  the  screen  downcomer;  however  the  main  bank  fluid 
absorbs  additional  energy  from  the  desuperheater  located  in 
that  circuit.  Therefore,  the  main  bank  riser  inlet  enthalpy 
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must  be  calculated  separately.   Thus, 

The  downcomer  density  must  be  calculated  for  use  in  the  pres- 
sure drop  calculations.   Hence, 

fcfkXO^   -   \..<f>  /\ft£xfi  \jrcf] 

The  riser  outlet  quality  is  calculated  using  the 
initial  quality  formulation  previously  developed.   This  com- 
putation is  followed  by  the  calculation  of  the  riser  non- 
boiling  and  boiling  lengths. 

The  average  density  of  the  risers  is  calculated  along 
with  the  two  phase  flow  multiplication  factors  for  use  in  the 
flow  rate/pressure  drop  calculation. 

An  updated  flow  rate  is  now  computed  and  compared 
with  the  first  approximation.   If  it  is  within  a  specified 
error  criteria  calculation  stops,   If  not  the  previous 
approximation  is  updated  and  the  calculations  continued  with 
the  new  approximation. 

Upon  completion  of  the  flow  rate  balancing  the  initial 
flow  rates  in  the  downcomer/riser  loops  are  known  along  with 
the  riser  outlet  quality. 

The  steady  state  drum  specific  volume  and  density 
are  computed  as 
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The  initial  steam  mass  in  the  drum  is  then  calculated  as 

9 .   Superheater  Pressure  Drop 

All  entrance,  head,  and  exit  losses  are  considered 
negligible  in  the  superheater  compared  to  the  frictional 
pressure  drop. 


AP     -       U 


A   ^average. 

By  lumping  the  constants  the  pressure  drop  equation  is 

Pmrv\-PwM  =  KckjI   ^M„  %  ^  ^  /RHcmiOjZS 

The  following  definitions  apply  here: 


10 .   Steam  Valve  Equation 

The  flow  through  the  steam  valve  is  considered 
directly  proportional  to  the  outlet  pressure  and  valve 
opening,  i.e., 
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rvwwxn:    -  pm<u  *  KorJ4    #  vavu/e,      \y\z\ 


The  constant  is  computed  in  the  initial  condition  program 


E.   EQUATIONS  OF  STATE 

The  equations  of  state  listed  below  are  used  in  both  the 
initial  condition  and  dynamic  programs.   With  the  exception 
of  the  subcooled  specific  volume  equation  used  for  the 
feedwater  entering  the  drum,  they  are  reasonably  accurate 
in  the  300-1500  psi  range.   The  subcooled  specific  volume 
equation  is  accurate  in  the  &0O  "*  %CC    psc  range. 

Ps/vr  -  exi° ((  Alos  (wskt)-  4  A<*70t)/. 2<q*4S2>         ]XJ 
BSAT  ^  HXP  (0.8fe4£2,  *  AL06  (psat)  +-  4.H6«^r)    W 

*:P5Ar^:^  2.0—  M.fc^E-^?**  3.0      Y&3 

RttOF    -  t3.?  -0.01TO1  ^T5AT^  1.1*32.  E.-0r 

*-  tsat  *  *  i.<f>  ~  (0.r7?(Qe-fe  *  ts/vt 

v/66  -  .^icoj6^h??-*.0^^^0H(o  )fe"re© 

-?.l4zE-u*-nsa  **3.0  -hi-H^sie-re 

5fc  T£>6  ^  ^  q.0  -  Ll4?E-\0?  #;TI3e  ** -T-0 
-h  ?.C2>4&-Z0*TBC  ^^G.0 
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IV.   RESULTS  AND  CONCLUSIONS 

A .   GENERAL 

A  listing  of  the  initial  condition  program  and  the  dynamic 
boiler  model  program  is  given  in  Appendices  A  and  B  respec- 
tively.  The  initial  condition  program  output  data  must  be 
properly  formatted  for  input  to  the  CSMP-III  dynamic  model. 
In  addition,  because  the  dynamic  model  utilizes  the  liquid  in 
the  steam  drum  as  a  saturation  state  point  from  which  all 
other  state  points  are  derived,  two  of  the  initial  conditions 
must  be  modified  to  eliminate  a  discontinuity  between  the 
steady  state  program  and  the  dynamic  model.   The  equations 
involved  and  an  explanation  are  given  below. 

DMPtSL^  -=  (vold*/v\  *  Rttaubtf) /z.tf  £rc[ 

D/vvdH-L  a  \tJTGR.\-(pMOtiU0)  DO MOH-l)  LmC 

As  stated  previously,  the  enthalpy  of  the  liquid  in  the 
drum  is  considered  saturation  enthalpy.   The  dynamic  model 
must  begin  with  the  rate  of  condensation  in  the  drum  (MCOND) 
as  close  to  zero  as  possible.  This  is  a  natural  steady  state 
posit  ion  -  drum  liquid  and  drum  steam  both  at  the  same  pressure 
and  temperature.   In  order  to  insure  this  to  be  the  case 
DMASLO  and/or  DMDHLO  must  be  modified  such  that  the  initial 
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drum  enthalpy  (DH)  very  closely  approximates  the  enthalpy  of 
saturated  liquid  corresponding  to  the  drum  steam  temperature 
and  pressure  (PMM  and  TMM).   This  is  easily  facilitated  by 
the  use  of  the  CALL  DEBUG  statement  in  the  dynamic  program 
With  this  procedure: 

1.  The  initial  condition  program  is  executed  and  the 
initial  conditions  formatted  for  CSMP-III  use. 

2.  The  dynamic  model  is  executed  for  only  a  short  run 
time,  i.e.  5  seconds. 

3.  Observing  the  DEBUG  output  from  the  model,  DMASLO 
and/or  DMDHLO  are  modified  such  that  PSAT  and  TSAT 
equal  PMM  and  TMM. 

4.  The  model  is  reexecuted  for  a  short  run  to  check, 
MCOND  should  be  very  small. 

5.  Since  DH  is  implicitly  related  to  PMM  and  TMM,  the 
procedure  may  have  to  be  repeated  a  few  times.   The 
objective  is  to  force  as  many  of  the  "DERIVATIVE" 
terms  in  the  DEBUG  output  to  comparatively  small 
figures  as  possible. 

The  model  as  used  operates  satisfactorily  with  DMASLO  = 
8002.2  and  DMDHLO  =  3.96163E  06.   This  forces  the  initial 
MCOND  term  to  .35949. 

Because  of  the  CSMP  function  DERIV  used  in  the  program, 
the  model  is  extremely  sensitive  to  integration  time  step. 
Numerical  differentiation  is  not  a  desirable  function  to 
perform  in  a  dynamic  model;  however  attempts  to  explicitly 
differentiate  the  equations  concerned  failed.   One  solution 
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to  the  problem  is  to  "smooth"  the  derivative  function  by 
averaging  it  over  several  timesteps.   The  model  performed 
satisfactorily  with  a  fixed-step  integration  procedure 
(Runge-Kutta) ,  an  integration  time  step  of  .04  seconds, 
and  averaging  the  derivative  over  sixteen  time  steps. 

B.   OPEN  LOOP  RESPONSE 

The  open  loop  response  to  a  ten  percent  increase  in 
throttle  opening  is  shown  in  Figures  2-7.   As  expected,  the 
response  of  the  main  bank  circulation  loop  and  screen  bank 
circulation  loop  is  different.   This  is  in  keeping  with  the 
different  purposes  of  those  two  loops,  steam  generation  and 
furnace  screening  respectively.   The  effect  of  the  valve 
opening  increase  is  barely  noticeable  in  the  screen  circula- 
tion loop. 

The  swell  effect  is  not  noticeable.   Further  conversa- 
tions with  Mr.  Paul  Weitzel  at  Babcock  and  Wilcox  indicate 
that  the  percentage  of  "carry_  under"  is  not  a  constant  as  it 
was  treated  in  this  program.   One  to  two  percent  is  a  good 
starting  estimate  for  steady  state;  however  during  a  transient 
"carry  under"  mass  flow  rate  is  computed  by  subtracting  the  amount 

of  steam  leaving  the  boiler  from  the  amount  of  steam  produced. 
This  should  always  be  a  positive  number.   The  vapor  that  does 
not  leave  the  boiler  is  "carried  under."   The  following  program 
will  implement  this  "carry  under": 
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Modify 

the  section  of  the  dynamic 

model  titled  "COMPUTE  THE 

DRUM  SPECIFIC  VOLUME"  to 

read: 

PROCEDURE  CRYUND=FILTR8(MMM) 
IF(TIME.GT.0.0)GO  TO  45 
CRYUND=PCU* (MLL+MFF ) 
GO  TO  46 

45  CRYUND=(MLL*XLL+MFF*XFF)-MMM 

46  RISE1=DELAY  (250 , RISTIM, CRYUND) 
RISE2=CRYUND 

IF(VALUE.GT. .51)GO  TO  55 
RISE=RISE2 
GO  TO  57 
55         RISE=RISE1 
57         CONTINUE 

ENDPROCEDURE 

DDMDV=((MFF+MLL-MBB)*VF+CRYUND*VV. . . 
+MBB-P*VBB-(MCC+MGG)*VF+MCOND.  .  . 
*VFG-RISE*VV) 

Follow  the  above  statements  with  the  unmodified 
equations  for  DMDV. 

The  previous  procedure  was  not  implemented  in  the  present 

dynamic  model.   The  quality  formulation  in  the  present  model 

is  apparently  too  simplified  and  this  causes  a  shortfall  in 

quality  at  the  outlet  of  the  screen  riser  bank.   As  a  result  of 

this  shortfall,  there  are  instances  during  transient  operation 

when  the  model  is  not  "producing"  as  much  steam  as  it  is  "using" 

Mr.  Weitzel  suggests  using  twenty  node  finite  difference 

approximation  for  quality.   This  could  be  done  using  an  equation 

for  quality  such  as: 


where 


X   -         9 * 

X  =  quality 

Z  =  distance  up  the  riser  (ft) 
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C\  =  heat  input  to  riser  (BTU/s-ft) 

o 
A  =  cross  sectional  area  of  riser  (ft  ) 

3 

pf  =  density  of  saturated  liquid  (lbm/ft  ) 

V    =  velocity  of  liquid  entering  (ft/s) 
h    =  latent  heat  of  vaporization  (BTU/lbm) 

Simultaneously  an  implicit  equation  for  continuity  could  be 
applied  to  produce  a  balanced  mass  flow  rate  without  assuming 
that  the  downcomer  mass  flow  rate  remains  constant. 

C.   CONCLUSIONS 

The  model  presented  is  not  a  finished  model.   Further 
research  is  needed  to  successfully  implement  the  shrink  and 
swell  theories  presented.   When  the  previously  mentioned 
difficulty  with  riser  outlet  quality  is  solved,  the  model 
should  approximate  the  dynamics  of  a  wide  range  of  D-type 
marine  boilers  depending  on  the  initial  conditions  supplied 
to  the  initial  condition  program. 

The  initial  condition  program  develops  the  initial  condi- 
tions necessary  for  the  detailed  boiler  model  with  the 
relatively  scant  data  available  from  the  boiler  technical 
manual  shown  in  Appendix  C  and  a  small  amount  of  data  from 
common  engineering  handbooks.   Because  of  the  relatively 
small  amount  of  data  required  the  initial  condition  program 
and/or  the  model  can  serve  as  a  basis  for  a  boiler  condition 
monitoring  system. 
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D.   SUGGESTIONS  FOR  FURTHER  RESEARCH 

It  is  evident  that  the  model  suffers  by  using  the  CSMP- 
III  function  DERIV.   In  order  to  facilitate  the  explicit 
solution  of  the  continuity  equation  for  the  riser  banks  some 
form  of  differentiation  is  required,  either  numerical  or 
explicit.   An  explicit  solution  is  much  more  desirable  from 
a  stability  standpoint.   A  more  detailed  investigation  with 
fewer  assumptions  might  produce  the  correct  form  for  the 
explicit  differentiation  equation.   It  should  be  noted, 
however,  that  when  using  explicit  differentiation  the  model 
must  average  the  derivative  over  at  least  two  time  steps  to 
avoid  the  creation  of  an  algebraic  loop.   This  is  easily 
done  with  the  same  PROCEDURE  format  used  to  average  the 
derivative  in  the  current  model. 

The  use  of  small  perturbation  techniques  to  linearize  the 
model  and  thus  allow  a  state-space  representation  for  multi- 
variable  control  development  and  analysis  should  be  undertaken 
There  is  a  computer  code  available  locally  to  facilitate  this 
for  CSMP  models;  however  it  does  not  accept  the  CSMP  DERIV 
function  and/or  the  averaging  procedure.   A  means  of  bypassing 
this  problem  would  result  in  a  general  D-type  boiler  model  - 
applicable  to  a  very  wide  range  of  boilers  currently  in  use  - 
which  could  be  linearized  about  specific  operating  points. 
These  individual  models  could  then  be  used  for  the  development 
and  testing  of  multivariable  control  systems. 

Using  locally  available  optimizing  computer  codes,  an 
optimal  D-type  boiler  design  could  be  attempted  with  regards 
to  boiler  geometry. 
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A  more  detailed  analytical  and  experimental  investigation 
of  current  D-type  marine  boilers  should  be  undertaken  locally. 
Little  research  has  been  done  in  this  area  and  that  which  has 
been  done  has  often  been  based  on  either  an  incorrect  physical 
model  of  a  marine  boiler  or  on  a  nuclear  steam  generating 
plant.   A  starting  point  could  be  the  data  analysis  of  common 
horizontal  and  vertical  steam  separators  followed  by  an  opti- 
mum design  for  these  element s„   The  improper  and/or  mainten- 
ance of  drum  internals  apparently  grossly  affects  boiler 
operation  and  water  level  stability  during  rapid  transients. 
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APPENDIX 


C8/3C/7S       2C.55.ll 
IL£:    CONSTANT    F0RTP4N       =>  I 


A 


NAVAL    POSTGRAOLATE    SCHCCl 


THIS    PROGRAM    CALCLLATES    the    CONSTANTS     AMD    INITIAL     CDNCITHNS 
FOR     A    US    NAVY    D-7YPE     BOILER    GIVEN    THh    INPUT    DATA    OBTAINED    F"CW 
ThE     ECIL£R    T£ChNICAL     '"ANlJALi     STEA'-«     TABLES,     AND    STAND  AfO    ENGINEER- 
ING   HANGBCOnS.        THE     INPUT     DATA    ARE     ENTERED     IN    MA -IE     LIST     FORMAT. 
NA^ELIST    DESIGNATIONS    AND     INPUT    NEMQNICS    FJLLJw 


INCONl: 


OPPNT 

TOT3TM 

SFHSTM 

OSHSTM 

PCRUM 

SHOT 

SHOP 

OSHCT 

DSHOP 

ECONIT 

ECONGT 

AI~TMP 

CILTNF 

AIRFLQ 

0  ILFLC 

XCSAI3 

DPAFT 

TGA3SC 

TGASSH 

TGASMB 

TGASEC 

TSCRN 

TSPHTR 

HRFVCL 

FHAS 


OPERATI i 
TOTAL  ST 
SUPERhtA 
DESUPERH 

CRUM  PSE 
SUPcRhEA 
SUP  E«  HE 
OES'JPERri 
DESUPERH 
ECTlTi  Ll 
ECC.NCN  11 
AIR  IMLc 
CIL  T  E"P 
AI RFt_  DW 
TIL  FlO* 
EX  C  £S  S  A 
DRAFT  LJ 
FLUE  GAS 
FLLE  GAS 
FLUE  GAS 
FLUE  uAS 
SCREEN  T 
SU°EC'  FEA 
HhAT  RcL 
FURNACE 
ABS'JR  BM 


G    PO 
EAM 

TER 
EATE 
S  SU  R 
TEH 
TE' 
EA'E 
EATE 
E3     F 
EP    F 
T    TE 
t  c  AT 
R  £~E 
RAT 
IK     ( 
SS    ( 
TEM 
"EM 
TEM 
TE'" 
U8E 
TE=> 
E*SE 
HEAT 
G    AR 


INT 

FLO 

STE 

o   s 
t  ( 

OUT 

nuT 

R      ) 

r  n 

EE) 
EEJ 
MPfc 
t)RE 
(  L 
E  ( 
»ER 
IN- 
PFR 
PER 
PER 
PER 
WAL 
TIB 
|.< 
Ad 
EA  ) 


(PERCENT)        * 
«(LB/HR) 

AV    FLOJ     (Le/HR) 
TEAM    FLOW     (LB/HR) 
PSIG) 

LET    TEMDEPAT'IOE    (CEG-F) 
LET    PRE  SSU<?E    (PSIG) 
UTLET    TE''PE0',Tijoc    (JEG-F) 
L'TLET    PRESSURE     I  PS  IG) 
INLET    TEMPERATURE     (OEG-F) 
TEVPE-ATURE    (DEG-F) 
TO    REGI STcRS 


>J)TLET 
RATURh 

AT  5J  = 
3/HRJ 

LB/HR) 
CENT  ) 
H2G) 

\TURF 

ATUPt 

ATURE 

ATURE 

L 

E 

B 


IZEG-f) 


1ER     INLET    I  JEG-F) 


LEAVING 

LEAV  IN'G 

LEAVING 

LEAV  IMG 
TE^PERATU»E 
WALL     TEMPERATURE 
U/FR/CU     FT    F'.JQ  MA 


SCREEN    (GEG-F) 
SUPERHEATER     ICEG-F) 
*A1N    BANK     (QEG-Ff 

ErJNJMlZER     (DEG-F) 
(DEj-F ) 

( DcG-F) 

E    VOLUME) 


SJRBTION     (K6TU/HR/SQ    FT     RACIANT     HEAT 


NOTE:    "PERCENT"    VALUES    ARE    TO    BE     LCGGEO    AS    DECIMALS. 
$ SCREEN 


DTuasc 

LAVSC 
NTUBSC 
RHASSC 
MASSSC 


AVERAGE    INSIDE    DIAMETER    IF    SC  ~  EEN    TUBES     (IN) 

AVERAGE    LENGTH    OF     SCREEN    TUBE     (FT) 

NUMBER     IF    SCREEN    TUBES 

RADIANT    HEA"     APSCRPING    AREA    IF    FURNACE    SCREEN    (SO    FT) 

TOTAL    WEIGHT    OF     SCREEN    TUBES     (LB) 


SSPHTR 


AREASH 
MASSSH 


JMNBAMK 


OTUBMg 

LAVM3 

NTU3M3 

Msswe 

AREAMB 


TOTAL    SUPERHEATER    HEAr     TRANSFER    AREA    (SO 
TOTAL     WEIGHT    3F    SUPERHEATER    TUbES 


FT  ) 


SECON 


DTLBEC 

NPASSE 
NTU9EC 
LTU8EC 
MASSEC 


AVERAGE     INSIDE    DIA^ETE0    TF    MAIN    BANK    TUBES 
AVERAGE    LENG"H    OF    M*IN    RANK    TUBE    ( IN i 
NUMBER    3G     -1AIN    BANK     TUBES 
TCTJL    WEiGtiT    hg    MAIN    R4NK    TUQES    (LB) 
TOTAL    HEAT    TRANSFER    AREA     (  SQ    FT) 


.INSIDE    DIAMETER    CF    ECONOM  I2ER  TUBES  ( IN  ) 

NUMBER    "IF       TUBE    PASSES 

NUMBER    IF  TUBES    °E°    PASS 

AVERAGE  LENGTH  OF  ECONOMIZER  TUBE  (FT) 

THTAL  HEIGnT  OF  ECONOMIZER  (LB) 


(  IN) 


CONOOO  10 
CGN00  0  20 
C0NCCC30 
CONOoCAU 
CQN0GG50 
CDNCCC60 
C0N00070 
C0N00030 
CONCCGSO 
CONOOIOO 
C0N0G1  10 
C0N0012O 
C0N00130 
CONCC140 
CONGO  150 
C0N0016J 
CDNCC17C 
C0N0013J 
C0NC0190 
CON  CO  2  CO 
C0NQ02  10 
CGNCC22C 
CON  0  02 JO 
CGNCC2+J 
C0NCC2D0 
C3N00260 
C0M00270 
C3NCQ230 
C0M0G2QJ 
C0NCC3C0 
CON  00  2  10 
CON00320 
CONCC330 
CONGO  _>  AG 
C0NG0350 
C0NCC360 
CONOG3  7G 
CONCC350 
CON  00 390 
CONOOIOO 
COM  CCA  10 
C0NC0A20 
CON0OA3O 
C0NCC44C 
C0N00450 
C0NG0^60 
CON  CCA  70 
C0N00480 
C0NC0A90 
C0N0U50U 
C0NCC5  10 
CONCC520 
CON005  3G 
C0N0G540 
CONCC550 
CON005oO 
C ON C 05 70 
CON  00 5 60 
CON005  90 
C0NCC6CC 
C0N00610 
CJN00620 
CONCG630 
CDNG06A0 
C0NCC65C 
C0MC0660 
C0N006  7G 
CDNCCfcEO 
C0NCG6SG 
C0NC0700 
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FILE:    CONSTANT    FORTRAN       PI  NAVAL    POSTGRADLATE  SCHCOl 

C  CONC071G 

:  CON007 20 

C  SCES°F  CGN00730 

:  concc74o 

c  dturds  inside   jiameterof   desuperhea  ter   tube    (in)  con0075c 

C  NTUBGS  N'JMBbF     IF    OESU  "ERhE  ATER    TUqES    PER    ^ASS  C0NCC76U 

:  NPASSO  NUMBER    CF    DESUPfcRHEATER'  TUBE    PASSES  CUNC077G 

C  LTII8CS  LENGTF     IF    3ESU=»ERHE  *TF3     TL3E    (FT)  C0NG073J 

;  APEAOS  T:TAL    HEAT    TRANSFER     A?EA    OF    OESUDERHE  ATE'     (SO  FT)                 C0NCC7SG 

:  MASSOS  TJTAL     HEIGHT    CF    Dt SLPERHE ATER    ASSEMBLY      (Le)  CON00300 

C  C0NCG8  10 

Z  CDNCC620 

C  iDRMOCR  C0N008  3J 

C  C0N0064G 

:  DTU3DD  AVERAGE    DIAMETER    CF    DRUM    CCWMCOMER   TUBES     (IN)  CON00850 

C  LAVCC  AVERAGE    LENGTH    OF    DRUM    OQWNCQMER    TUBE     (FT)  C0MCC860 

C  NTL8CC  NUMeER    OF    DRUM    CCWNCOMER    TueES  CONCC67Q 

C  -  C0NG0830 

C  C0NCG89U 

C  SHORDCP  CON00900 

C  CCN0G910 

C  DTUBHC  AVE' AGE    INSIDE    DIMETER    OF    HEADER    DOWJCOMER    (IN)  CON00920 

C  LAVHD  AVERAGE    LE^TH    CF    HEAGER     C^-J  NCOM  ER     (FT)  CONGC930 

C  NTU3HD  NUMBER    JF    HEADER    00  WC  3  ME  R*S  CGNQ0940 

C  CONCC950 

C  CCNOG960 

C  saCILER  C0NGC970 

Z  CONCG9(3G 

C  OST^DM  DIAMETER    OF    STEAM   DRUM    (IN)  ■                         CON00990 

C  LSTMDM  LENGTF    IF    STEAM    qpum    (FT)  CQN01000 

G  DWTRDN  DIAMETER    CF    WATcf-     D  F'J M     (IN)  C0NC1010 

C  LWTRDH  LENGTH    JF     WATER    DRUM    (FT)  C0N01020 

C  HN'Ci^M  hcluHT     ^F    NCRTAL    WATER    LEVEL     ABOVE    BENCH    VARK  (FT)                CONC1020 

3  HHDR  HEIGHT    CF    HEADER     (SCREEN)     AGGVE    BENCH    MARK     (FT)  CON01G4Q 

C  hwTRCM  FEIGhT    IF    »'.*T£_*    DRUM    ABOVE    BENCH    HARK    (FT)  CGNGU5Q 

C  FCRVCL  FLPNACE    VOLUME    (CU     FT)  CONCLC60 

C  CO  NO  10  70 

C  NCTE:  CHOICE    OF    3ENCH    MARK. IS    ARBITRARY  COMOL080 

C  C3NCIC90 

C  JTHERMO  ..  CONOilGJ 

C  CONGlllO 

C  HJHOLT  ECONOMIZER    CUTLET    ENTHAL0Y     (BTU/L3M)  CONC1120 

C  HDSOUT  OESUPtRHEATER    OUTLET    ENTHALPY    (3TU/LBM)  CQN0113J 

C  HECIN  ENThALPY     CF    ECTn^lZE^    FEED     INLET    (3TU/LBM)  C0NC114J 

C  HECOUI  EiJTHaLPY    OF    ECONOMIZER    FEED    OUTLET     (eTU/LEM)  C0MC1150 

C  (DKh2:  ThE^MuL    CONDUCTIVITY    JF     *ATER        (  6  TU  /HR.  F  T.DEG-F )  C0N0116G 

C  IDFPH2C  PRANGTL    NUMBER     FOR    W  AT  E  R  C0NCU7G 

C  (DVSrH20  KINEMATIC     VISCOSITY    FOR     WATER     (  LB  VFT-SEC  )  C0N0U80 

C  (2)KST*  THERM.iL    CONDUCT  IV  ITY  FOR     STEAM    (  3TU/HR  .FT  .DEG-F  )  CGNO  1  190 

C  (2JPRSTM  PRANDTL    NLM3ER    FOR    STEAM  CON01200 

C  (2JVSC3TM  KINEMATIC    VISCOSITY    FOR    STEAM    (LBM/FT-SEC)  CGNC1210 

C  R5HCLT  SUPERhtAT£R    CUT  LET     GENS  ITY     (LE^/C.J     FT)  C0N0122O 

C  RDSOUT  0E5UPERHE ATER     ULTLET    DENSITY     (L3V/CU    FT)  C0N01230 

C  RFLUE  CENSITY    OF    FLJE    GAS    AT    TGASSC    (LBM/CU    FT)  C0NC124C 

C  CONC125G 

C  NOTES:  (1)     EVALUATED    A7    AVERAGE     ECONOMIZER    W^TER    TEMPERATUCE       C0NG1260 

C  12)     EVALUATED    AT     AVERAGE    DESU^ERHE  ATER     STEAM    TE MPERA TUCCUN C 12 7G 

C  C0M01230 

C  $LCSSES  C0NG129J 

C  KSCSC  ROUGHNESS    PATION    (SAND    EQUIVALENT)     FOR    SCREEN  TUBFS            C3NC13CG 

C  KSDMB  ROLGHMESS    RATIO    (SAND    EQUIVALENT)     OF    MAIN    BANK    "HJBE            COM01310 

C  KSCCC  ROUGHNESS    RATIO     (SAND    EOU  IV-L  ENT  J    FOR    DRLM    DOWNCOMEP         CGNC1320 

C  KSDHO  ROUGHNESS     RATIO    (  Sa  fO    EQUIVALENT)     FCR     hEACER    DO^NCCER    CON01330 

C  ENTSC  SCREEN    TU3E    ENTRANCE    LOSS  CONCUR 

C  EENCSC  SCREEN    TUBE    BEk"0    L~SS  C3NC1350 

C  ENTMB  MAIN     BA  JK    ENTRANCE     LCSS  C0N01360 

Q  BENCMe  MAIN     feriK    BEND    LOSS    FACTOR  CON01370 

C  EfaGG  DRUM    CCWNCCttER    ENTRANCE    LCSS     FACTOR  C0N013S0 

C-  BENDOO  DRUM    OC.'JNCGMER    3E  ND    LOSS    FACTOR  CON01390 

C  ENTHC  HEACER    CCWNCOMER    ENTRANCE   L">SS    FACTOR  C0NC1AC0 


09 


CONST  *MT    FORTRAN       PI 


NAVAL    POSTGRADUATE     5CHCCL 


8ENDH0      HEADER    OOWNCCMES    BENO   L^SS    FACTOR 

IMPLICiT    REUIA-Z) 

INTEGER     ISPASSE 

DIMENSION    TH2J(20),OOASS(2C) 


OEFINE    N4MELISTS 

NAfELlST/INCCNO/CP 

lDShT0,EC3  II  T,EC3N  ) 

ZTC-JSSF.-CASVR.TGAS 

NAVELI ST/SCRtEM/DT 

NAMEL  IST/SPHT'/V^F 

hA*»EUST/VNE  dNK/  LT 

NAMELI  ST/ECCN/JTLB 

NA"EL  IST/CES^H/DT'J 

f^Af-ELI  ST/Gh"DC^/LT 

NAMEL IST/HORGCR/DT 

NAfrELIST/BOILER/QS 

SFLRVCL 

NAMELIST/THERM1/HS 
lKH2CtPRH2C,V^  C 

NAMEL I iT/L 1SSES/KS 
iBENC^B,  E  J7C0,  -3E  i  CC 
MPELIST/INCCMl  /TR 
6TYYCTZZC,  I  AXJ,7^ 
NAMEL  IST/  INCCN^/f-'A 
SMNNOfMMMO.MQQO t 
£MCDO,MGH0,  •nr.OO, 
GMFFO  ,NLLJ  ,MVE, 
€HASSV,l«A  SSU  ,,vASSY, 
EMASSXV  -lASSZt^ASSOS 
EO^UHIC, JNASL3, JSTv 

tLSTMOM  ,  J  ST  1D.M 
NAyELI3T/ INCONi/ f*v 

ccpy ,cpx  ,cpz» 
ecpwtco^  ,crs,:st, 

SCTU  ,Cf>F,  J,GC 

NAMELI  ST/CONST1/KS 
6KSCC»KS0E»KSEFf 

SKSJK  ,kskl ,KZ  t 
£KW, <0S,KS  r, 
CKQM  ,KTIJ,  KCN4, 
£KX,KHJ, K\, 
£KTU,KY 

MPELIST/C0NST2/  £  J 
SLJLltDJF.LDF  , 
CLCFl,  ACF,LGF, 
£AGH  ,CGH»CCG  t 
£LCG,LD£G,  LJ<<G 

NAWELIS"/C0NST3/  jc 
6  VOLJL  ,VCDF  ,VCLOhN 

£VOLHJ,ZOF,Z JL, 

S.LJL  ,DDF  f.jKLt 

6CEF,0JK,ZCD 

NAVELIST/CONSTV  FC 
6FEF  ,FJis  tFKL. 
*ENTRGH,E  ITR  COf  BEND 
6  8ENCCOfEXl'JH,EX  IT 

NAMELI  Sf/C  nST3/FH 
£TA"B,C5TCFC,  DROJLCI 
£SIGf*£ 

NAMELIS7/OU1PLT/Q 1 
£XFFG  ,XL  LOt*  FCGFO,  P 
6HSATtHCDC  ,HAAG  tHEM'. 


PN  T,  TOT  S  TV  t  S  FH3  T  M,  CS  HS  T'«,  PORJM,  SHOT,  SHIP,  DSFrT, 
rtMRTMP  ,JILFLOiAIRFLC|XCSAI  RtURAFT  ,TGASSCi 
EC  fTSCRNtTSPHT^f  HP  FVTLt  FHA  i,  OIL  TMP 

03  SC  tLAVSCtNTUBSCf  PHA5SC,M  AS3S  C 

\SH,  -'ASSSH 

J  B  v  B,  L  AV  M3 ,  M  TU9Mof  m  ASSM  B ,  AS  £  AM  R 

£C  ,NPASSE,NTUbEC  ,  LT»  'BE  C  ,  MAS  *  EC 

^3,  MTLPt)S,•|P,.  SST.LTUSOS  fARE  AO  S  ,  »  A  S  SO  S 

!JeCCtLAVCDt\TUBOD 

UGHO  ,LA VHG  ,N7U3HD 

"  M  C  t* ,  LS  T  '-'J*,  OWTP  OM  ,  L  WT  ?  DM  ,  r-HOR  M  ,  HHDR  ,  HWTRD  M  , 


<ST 
OS 


CCJ,  FN  T 
TRR':  ,T\ 


lT,HDSOUT,HECIN  ,HECOUT, 
M,  PRSTMfV5CSTM,  RSHOUT,RDSO!JTtRFL 
,K^J MB  ,  is SOOG  ,  KS  GHJ  ,ENTS  C  »  6 EN CS  C  , 

MTHD,RENDHO 


UE 
ENTMB, 


wo  ,twwo 


iO»MFOQO,MAOQO« 

)  ,C$TMDM, 
Nt CAB.CPV, 

JL,KSOF  ,KSGH  , 

L,CJL,LJL2, 

CtZBENCltZBENC2, 

» 

0,  FGH,  FOE, 

jHf 

CC 

V, >LL0,XFF0, 


G  ,C40,Q6C,QeC,Q90,*FFO,MLLO,  CCO  , 

hOJLJt  RH-tFFOt*  hOLLO, 

J  fhJJO  ,HFGfRhCFf  OMCVO,  HXFFG,  hXLLG 


Ac 


Pc^0(5  ,  INCCNO  ) 
REAO(  5,  SCREEN) 
PE4C(5,S°hT3 ) 
READ  (5,  1^8  ANK) 


CONOID  10 
C0N01420 
C0NC14  3C 
C0NG14A0 
CON01<t5G 
CON  C146C 
C0NG1470 
CON014  80 
CJMC149G 
C0N01500 
CGN0L510 
C0N01520 
CCNG15GG 
CONCIS^O 
CONG155G 
C0N0I560 
C0NC157G 
C0NC158G 
C0NC1590 
CONC16C0 
C0NG16LG 
COM C 1620 
C3N01630 
C0NC1&40 
CON  C  it  50 
CONG  16 60 
C0N0U7C 
CONG  16  80 
CCN016SG 
C0NC17CG 
C0NC1710 
CO  NO  17 20 
CONC173C 
C0N0174G 
COMC17  50 
C0NC176C 
C0N017  7G 
C0NC1 76C 
CONG  1790 
C0NC1800 
CONOiaiQ 
C0M013  20 
C0N01c3C 
CONG  18  AO 
CQNC185G 
C0NQ1860 
CON0187J 
C0NC183G 

conci8<;o 

C0N019G0 
CON019  10 
C0NC1S2G 
CONG  19 3G 
CON  019  AG 
CQN0195G 
CON C  19  60 
CONC197G 
C0N0198G 
CONC1990 
G3N0  2000 
C0NC2010 
C3NC2020 
C0N02030 
C0N0  2CAC 
C0NC205O 
CCN02060 
CONC2Q7G 
COM C2 0  80 
CON02U90 
C0NC21C0 
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FILE:    CONSTANT    FORTRAN      pl 

R£AD(  5fECCN) 
BE  JSC  (5,  CE5nM 
REAG  (5  ,CM:jCaj 
RE">C<  5,H^.)C?  ) 
PEAO  (5,  6nlL  f  1 
.     READ  (  5  ,THckVG1  • 
3EA0< 5, LOSSES) 
C  CALCULATION    CF    PRELIMINARY     CONSTANTS 


NAVAL    POSTGRADUATE    SCHOOL 


FI  =  3. 1415*27 
G=32.2 

PR AA=P3H23 

VISCCA=VSCI-2T 

1HC0NA  =  KH2  J 

pq  AN  =  PR  ST^ 

VI  <=CCN  =  VSC">T* 

THCC.NN  =  KS7^ 

XASUME=O.G5 

FC  L  =  .0I5 

GC  =  22.2 

VALVEO=3PP-JT 

KSCD=KSOhO 

KSDF=KSJSC 

KSDEsKiJF 

KSEF=KS3F 

CM*. 11 

CPV=0 

CP  Y=CM 

CPX=CM 

CPZ=CM 

CPW=CM 

KSGH=KSGCO 

KSJL  =  K50NB 

KSJK=KSJL 

KSK  L  =  KSJL 

ENTRCD=E^rnD 

ENTFGh=  ENTDC 

9EN0CD=*cNGH:) 

3ENDGF=3ENJ  JH3 

EXITCD  =  l  .0 

EXI  TGH  =  I.  J 

PS  AT=PC4UM+l4.7 

HSAT=EXP(0.2&45P*aiCC(PSA~)*4.46  703) 

HFG»S22.15-C.4C516*PSAT+1.  7  I  7E-04*  PS  AT**2.  0-4.  219£-0  9*PS  AT**3  .0 

FV*hSAT+hFG 

HMMO=H  V 

FF=hStT 

TS  AT=  EX  P  (.  22  15  1^  A  L  Co  (OS  4T  j  +4.7  7123  ) 

MFCQC=OI  LF|_G/3oOl.  J 

MACCO=  t  HFL"i/26oo.J 

7RF0=ToASSG 

TVVG=TSCVJ 

TV.WO=TSPhTR 

TQRO=TKRC 

TSSO^CNSSH 

THK)=TS  AT 

TTTO=TGAS^B 

TUUO=TGASEC 

TTLG=(TT10+TUU0) /2  .0 

TA*<3=8C.C 

*MM0=T  3TST-1/3600  .0 

MMNE=M*MC*.000l 

VMMHI=SPHST>'/J60G.0 

TNNO=SriCT 

TPPC=DSHOT 

PNN0  =  SF")P*14.7 

F»»N)*FDRLH 

MMM  I  I»JShSTM/360C.a 

MNNO=MV^I I 

TAAO=ECGM  T 


C0N021  10 

C0N02120 

C  TIN  C*  130 

C3NG2140 

C0NC<:153 

C0NC2160 

C0N02170 

C3NC21 tG 

CON G2 I  90 

C0NG22  0G 

C0MC22  1Q 

C0NU2220 

CON022JO 

CON  G2240 

C0N02250 

C3NC2260 

C0NC22  70 

CON02280 

CONC2290 

CCN02300 

CO  MO  2  2  10 

C3NC222G 

C3N023  30 

C0NC234C 

CON  C2 3  50 

C0NC2360 

CONC23  70 

C0N02360 

CO  NO  2  390 

C3N024Q0 

C0N02410 

C3.NC24  20 

C0NC243G 

C0N02440 

C0NC245G 

C0NQ2460 

CON02470 

CGNC2480 

C0N0249O 

CONC25CG 

C0NG2510 

C0NC2520 

CONC2530 

CCN0254J 

CON025  50 

CON02560 

C0N02570 

CONC2560 

C0NC2590 

CON02600 

C0NC26  10 

C0N02620 

CONC2620 

C0NC264G 

C0N02650 

CnNC266G 

C0NC2570 

C0NG2680 

CONC26SO 

C0NJ2700 

CONG271G 

CON02720 

CON02730 

CON C2 740 

CON027  50 

CON02760 

CON  027 70 

CON02730 

CONC27SC 

C0NG2SCQ 
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FILE:    CCI^STAfvT    FCPTP-AN       =>1 


NAVAL    POSTGRADUATE    SCHOOL 


C 

c 
c 

c.s 

c 
c 


TRR0  =  EC3vlJT 

VBeQ=.01600488-.000a320l46*T3  SO.  000  0000  36  51  I*  T^R  J**  2.  0 
G- 8.  1 42c-  lL*T^tt0**3.0  +  l.4O>*  L  E-L  3*T  6  30**4  .0-1  .14  9  E- 16  * 
fiTBRO  **b  .0*o.u  3  4E-  <J*T«?u**6.C 

CTCT0=n3F7CL*L00  0.0/jo00.0 
.    HNN0=H5HGUT 

CPJR*e=FFAS»lOOO.Q/3tOO.O 

AR4C=RH£SSC 

HPPC=HUS1LT 

LCF1=LAVSC 

L0F2=L3cl*NTUBSC 

CDF=  0TJ3  SC  *M TUB  SC/ 12.0 

C0E=CCF 

0EF=0DE 

A C F=  P  1* C7< I  3SC** 2. C*NTinSC  /(  4.0*144.  0) 

VCIDF=ACF*L0F1 

MASSV=MiSS3C 

WASSV=*A5S3  H 

LJL1  sLAVI"8 

LJL2=NTUB^3*LJL1 

OJL=DTiJa*a*NTUdVfc/LZ.O 

DJK=DJL 

DKL'DJK 

AJL=PI*CT'ie^E**2  .0*NTUB'W  (4  .3  <144.0  )• 

V0LJL=AJl*l_JLl 

MASSY=^3SS'13 

►ASS>*MA$Sfc:C 

MASSZ=1A,5i0S 

ZBENC1»M1RM-HHDP 

LDF=ZEtNCl 

ZDF»Z«ENC1 

ZeEMC2  =  FNCRf-hWTR0M 

LJL=ZBENC2 

ZJL=ZBENC2 

ZCC=ZCF 

ZGH=ZJL 

CR0CF0  =  0  .0 

0RCJL0=0. J 

LAB^LTUBEC 

CAe=cTueEC/  u.o 

AREAEC  =  NTU:3EC*LAc*ei*DAB*NFASS  E  "" 

DNPsCTUBCS/  12.0 

lCD=LAVh£*NTU  eHO 

LGH»LA.  \>00*NTU8D0 

DCC*0TUSFQ**ITU8HC/I2.a 

D(iH=L;Tuecc-'jruaoc/,i2  .0 

ACO*NTU3H3*(?I*{  0TURH0/24..  0)**2.0J 

AGh="ITUeCC*(0  I*(  Cr:i°00/24.u)  **2.0) 

V0L0FM=(  PI*  J STMOM** 2.0/  (4.0*144.01  )  *LST"CM 

VOLHJ=(0  [*0WTR0M**2. J/<  4.0*144) ^LWTRCM 

HAflO=HEC  IN 

H360=HECCLT 

RHOF=63.d-O.0173  1*TSAT+1.122E-C5*TSA  T**2.0-6< 

VF=1  .0/ fihCF 

\.FG=S  24.  C/PSAT-0.  1 

W=V  F-*VFG 

RH0V=1 ,0/VV 

RHOFG=RH3F-:'H'JV 

PHQM*0  =  3hCV 

RHONN0=R  SH3LT 

RH0PP0=RCSliJT 

MAA3=N"MC 

M8SC=MAAC 

CALCLLATE    TbE    TCTAL    MASS     CF    FLUE    GAS     IN    ThE    FURNACE 

MASSCR=RFLUE*F>JRV3  L 

CALCULATE    THE    TOTAL    ENERGY  ENTERING    BOILER 


7S6E-03*TSAT**3.0 


CCN02810 

CONC28  20 

C3NC2fc20 

CQN02840 

CONC2850 

C0NC2860 

CON02870 

C3NC268G 

CON02890 

C0NG2900 

C0N029  10 

CON02920 

C0NC2930 

C0MC2940 

CON02950 

C3NC2960 

C0N02970 

COMC2980 

C0NC2SS0 

C0N03000 

C0NC3CL0 

COM03O2O 

C0NC303O 

CONC2C40 

CON03O5O 

CON03060 

C3N03G70 

C0N030  80 

C3NC3090 

COM  03  ICO 

C0N031  10 

CONC312C 

CON  03 130 

C0N0314O 

CONC3150 

CON03160 

C0N0317J 

C3N0313J 

C0N03190 

CON03  200 

C3NC3210 

CQM032  2J 

C0NC3230 

C3N03240 

C0NC325J 

CON  03260 

C0N0327O 

CO NO 32 30 

C3N03290 

C0NO33OO 

CONC3310 

CONG  33 20 

CCN03330 

C0NC334C 

CON  033  50 

C0NC3360 

C3NC3370 

CON03380 

C0NC33S0 

CON  03400 

CQNC34L0 

C3NC3420 

C3N03430 

C3N03t40 

C3NC3450 

C0N034&0 

C0NC34  70 

C3N03480 

C0N03490 

C0NC35C0 
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NAVAL    POSTGRADUATE    SCHC01 


FILE:    CONSTANT    FORTRAN      »1 
OQOCTOTC*FURVCL 

CALCULATE    ••UI'IPEC"    LOWER    HEATING    VALUE    OF    FUEL/AIR    MIX 

FHV=CCG./*FOGG   ■ 

CALCULATE    ENERGY    ABSORBED    BY    RAOIATICN    IN    FURN4CE-C 
C1C=CFIF*B*ARAC 
C2  0=C10 


C 

c 
c 
c 

c 

c 
c 


c 
c 
c 
c 

c 
c 
c 

c 
c 
c 
c 


c 
c 
c 
c 

c 
c 
c 


c 
c 
c 
c 

c 
c 
c 
c 


CALCLLATE    THE    STEFAN-BCLTZ  *AN    CCNSTANT    MULTIPLIEC    BY    THE 
SCREEN    AREA 

SIGMAA  =  Q  10/(  (  TQRC+46C.0)**  A.  0-(TVV0+460.  C)**4.0) 

CALCULATE    THE    TCTAL    FLUE    GAS    '<ASS    FLOW    RATE     INTO    TOILER 

fCCO  =  NFCCJ*w/5CCO 

MRRU=MQ^C 

VSSO^'^C 

MUUC=*TK 

CCMFUTE    THE     ENEPGY    TRANSFERRED    TO    THE    W AT E3     (LOWER)    DRUM 
BY    THE    OESUPERHEATER 

Q90=NNNG*(HNN0-H PPO) 
G99G=C9J 

CALCULATE    THE    ENERGY    TRANSFERRED   TO   the    MAIN    SANK    RISERS 

C6C=K^O*  (HMMO-hEE3)*C20-0990 
05C=O6G 

CGMPLTE    the    VAfN    BANK    HEAT    f RANSFE3    COEFFICIENT-TUBE    yETAL 
TO    STEAM    SIUE3 

KY=AREAM3/1.  782EC6 

CCVF'JTE    THE    **AI\     BANK    TUBE    METAL    TEMPERATURE 

TYY0*(<:d0/(KY*PSAT**(<i.C/3  .0  )  )  )**(  1  .  0 /3.  0)  ♦■  T5AT 

COMPLTE    THE    MAIN    BANK    HEAT    TRANSFER    COEFFICIENT-FLUE    GAS 
TO    TUBE    M£T£L 

TSTO=(  TSS0+TTTGJ/2.C 

KST  =  C5  0/  ('•'SSG**0.o*(TSTG-TYYG  )  J 

COMPUTE    THE    SCREEN    R-4NK    HEAT    TRANSFER    COEFFICIENT-TUBE 
ftlAL    TO    STEAM    SICES 

KV  =  C10/  (PSAT**(<t.G/3  .0)*(  TWO-TS\  T)  **3.  0) 

COMPLTE    THE    HEAT    TRANSFER    TO    THE    SUPERHEATER 

04C=MMMG*(HNND-HV^J  ) 
030=040 


CCMPUTE     THE     SPECIFIC    HEAT    CF    THE 
IN    TJ-E    SU=>E'  HEATER 

CMN=OAG/(MMMG*(  TNNC-TM^O)) 


;team 


COMPUTE    THE     «;UPEPHtATER    HEAT    TRANSFER    COEFF  IC  I ENT-FLJ  E 
GAS    TO    TUBE    '^TAL-C 

TRSC=(TRRC+TSSO)  /2.  C 


CON  035  1G 
C0NU352J 
C0NC3  53C 
C3N03540 
C0NC3550 
CON C 3560 
CON  03 5 70 
C0NGJ58J 
CONC3  5SG 
C0N03oUG 
C0NG3610 
C0NC362G 
C0NG3630 
C0N0364G 
CONG365G 
C0N03660 
C0NC26  7G 
CONC366G 
CONG  36^0 
CONC37CC 
CON G3 7  10 
CON03720 
COTJ  C37  30 
C0N0374U 
C0NQ3750 
C0NC3760 
C0KG377G 
CONC27  6Q 
C0NC375G 
C0NG3800 
CONC381G 
C0N03820 
CONC3830 
CONC3640 
C0N03S5G 
C0NG3660 
CONQ3870 
C0\03880 
CON C 28  SO 
C3NC39C0 
CONC3910 
CONC2S2C 
C0N03930 
C0NC3  9tG 
C0NC3950 
CCN0396G 
CON03970 
C0M0398G 
C0N0399G 
CONC4C0G 
C3NC4010 
C0N04G23 
CONC4C30 
CON  04040 
C0NC4G50 
CONO^tGoO 
CONG4070 
CONCACcG 
CON  040 9G 
C0N04100 
CONC4U0 
C0NC4120 
C0hQ4L30 
C0NC4  14C 
C0NG4  150 
C0NC416C 
CONG4  17G 
CQN0418G 
C0NC4LCC 
C0N04200 
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FILE:    CONSTANT    FCRTRAN       PI  NAVAL    POSTGRACbATE    SCHOOL 

KRS=030/  (MRR0**O.  6*1  TRSO-T^WOJ  )  C0N04210 

C  C0NC422C 

C                            COMPUTE    THE    HEAT    TRANSFER    TO    THE    ECCNOMIZER  CONG42JG 

C  C0NQt240 

.     C8C=N  A  AO*  {he-30-h/541  J  C3NC425G 

Q7C=08G  C0N04260 

C  C3NG4270 

C                            COMPUTE    THE    SPECIFIC    HEAT     CF    THE    FEEDWATER  COMG42dG 

C                             IN     THE    ECONOMIZER  C0NG4290 

C  CONC43C0 

CAB=Q8C/( «AAC*( TBBO-TAAO) )  C0NC43 10 

C  CON04320 

C                            CCffLTE    THE    SPECITIC    HEAT    qf    FLUE    GAS    IN    THE  FURNACE                        CQNC433Q 

C  CQN04340 

CQPMCQ0-Q10  )/  (MR^O*(TQ*0-TAMB  ))  CON C4  35  0 

C  C0NC4360 

C                           COMPUTE    THE    SPECIFIC    HEAT   OF    F L IE    GAS    IN    THE  SUPERHEATER              CCN04370 

C  CONC43  6G 

CRS=Q3G/(MRRC*(TRRC-TSS0) >  C0NG43  90 

C  C3NG4400 

C                             COMPUTE    Th6    SPECIFIC    HEAT    OF    FLUE    GAS     IN    THE  *AIN    BANK                     CON0441J 

C  C0N04420 

CST=Q50/ C"SSG*(  TSSJ-TTTO  ))                                •  C0NC443C 

C  C IN  0  44 AG 

C                             COMPUTE    THE    SPECIFIC    HcAT    OF    FLUE    GAS    IN    THE  ECCNOMIZER                 CCN04450 

C  CONC4460 

CTU=07G/(  l»TT0*(TTT3-TUU0J  )  C0NC44  7G 

C  C0NG448Q 

C                            CCVFUTE    THE    ECCNC'UZEP    HEAT    TRANSFER    COEFFICIENT-TUBE  C3NC44S0 

C                            METAL    TO    LICUIu  C0N045C0 

C  C0NC45L0 

KX=(  U02  3*TrC1NAl/CAS)*(4.0/  (  P  1*0  A  3  *V  IS  C3A*NT'J  BEC  )  J**C.8  C0NC45  20 

6*PRAA**C.4*A:,EAFC  C0NJ453J 

C  C0NC454O 

C                             CCMPLTE    THt    LMTC    FOR    T HE    ECONOMIZES  CON04b50 

C  C0N045SG 

LMTDAe  =  08  3/  (KX*MiAG**G.3)  CGNC4  5  70 

C  C0NC45E0 

C                             COMPUTE    THE    ECONOMIZER    TUBE    METAL    TEMPERATURE  C0NCh!>9J 

C  C0NC46CC 

EXP0EC  =  E*P{  (  TBSG-TUC) /LMTDAB)  C0N0A61G 

TXXG=ITAJO-TeBo*EX30  tC  )/(  1  .o-EXPOG:  )                            •  C0NC4620 

C  CON C 46 30 

C                             COMPLTE    THE     ECONOMIZER    HEAT    TRANSFER    COtF FI CI  EN T-FLUE  CGNG4640 

C                             GAS    TQ    'USE    "5TAL  C0NC465C 

C  C0NQ46fc0 

KTU=0  7G/(MTTJ*(  TTUD-TaXO))  CON046  70 

C  C0NC4660 

C                            COMPUTE    THE    DESLPERHEATER    HEAT    TRANSFER    CCEFF  I  C  I  ENT -f'J  B  E  C3NG46S0 

C                           METAL    TO    S7EiM  C0N04700 

C  C0N047  10 

KZ=( (  .Q23*THCONN)/DNP)*< 4. 0/ < PI*DNP* VI SC0N*NTU3DS)  )  C0N047  20 

£**0.8*OP4N**J  .  i*ARE4.CS  C0N0  4  730 

C  C0NC474Q 

C                            COMPUTE    THE    DESUPEP  HEATER    L^  TO  CGN0475J 

C  C0NC476C 

LMT0NP=Q(;G/{KZ*MNN:**C.3)  C0N04  7  70 

C  CON04780 

C                            CCMFUTE    ThE    SPECIFIC    HEAT    OF   STEAM    IN   THE    OESU° ERHEATEP  C0NC47S0 

C  C0N048  00 

CNP=C90/  (MNNOMTNNJ-TPPO  )  )  C0NC4810 

C  C0N04320 

C                             COMPUTE    DESLPERHEATER    TUBE    METAL    TEMPERATURE  C0N0483G 

C  C0NC484G 

EXPOD$=E>P(  ( TNNG-TPPC)  /LMTCNP)  C0NC48  5C 

TZZO=  (TNNJ-EXPOOS^TOPO/I  1  .G-E  XPDDS)  C0NC4860 

C0NC4£7C 

COM°UTE    THE    DESLPERHEATER'  HEAT    TRANSFER    CCEFF  IC  I ENT  -  C0N0488G 

C                           WATER    OPUM    LIOUIO    TD    TUBE    METiL  C0NC4890 

C  CQNC49C0 


8 
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fILE:     CONSTANT     FCP-TRAN       PI  NAVAL    POSTGRAOLATE    SCHOOL 

KHJ=CS9C/(  TZZC-TSAT)  C0N04910 

C  CON0492J 

C  CGMFLTE    The     SUPERHEATER    L»-T0  CONC493Q 

C  C3N04940 

LMTCF»N^(TNN'0-T'4v|G  )/  (  *L1G  (  [TWHO-T^MQ  J  /(  TWViO-TNN  C)  )  )  CONC4S5C 

C  CON  0^9  60 

C  CALCULATE    THE    SUPERHEATER    HEAT    TRANSFER    COEFFICIENT-TUBE  C0NC4970 

C  MET*L    TO    STEAM  C0NC4980 

C  CONC4SS0 

KW=CKO/ ( ^MM j**o. a*LwTDMN )  CON05000 

C  C3NC5C1C 

C  CALCULATE    THE    THPJTTLE    VALVE    FLCW    COEFFICIENT  COMC5-020 

C  CON05-G30 

KGN4-MM«IIi/(\/ALVE)*F,\N0  )  CON  CSC  40 

C  C0M05O5J 

C  CALCULATE    THE    SUPERHEATER    OUTLET    DENSITY  CON05060 

C  CON05070 

•    K0N3=(?M'*0+PNNJ) /(  IHOMMO  +  RHONNC)  CON0508J 

PHCf  NC=  (PHG«yO+SI-CNNO  Ml  .0  C0NC5C90 

KQM-(  (  PMC-FMO^HCNhO)  /N'MMO**2.0  COM  C  51  CO 

C  CON0t.UO 

C  CCNFUTE    THE    F?ICTPN    FACTCRS     FOR   THE    OOnNCOMERS  CONC512C 

C  AND    THE    RISERS  C3NQ5130 

C  CONC5140 

FCD*l. 0/(1. 74-2.  0*ALCGlO(KlSCO)  J  •  CONC5150 

FCH=  1.0/ (  1.74-2.  C*\L0G10(<SGH)  )  C0N05-160 

FOE*  1.0/(1. 74-2. Q*ALCG10  (KSCE)  )  CONC517C 

FEF*1.0/(  I.  74-2.  C*AL3vilO(KSfcF)  )  CON05-130 

FJK=1.0/  (1.74-2.U*M.  1GU(  <\SJK)  )  CON05190 

FKL*  1.0/(1  ,74-2.G*ALCG10  (KSKL)  J  CON  C5  2  CO 

C  C0NC522C 

C  S7AH1    ITERATICN    TO    BALANCE    CIRCULATION    L01PS  CONC523G 

C  •  CONO!)240 

C  CnN0b26O 

C  COMPUTE     INITIAL    VALUE    OF     RISER    OLTLET    MASS    FLCVi    RATE  CON05270 

C  CON  C5  280 

HFFC=0  10/(  XflSUME*HFG  )  CON05290 

PLLO=C50/  (X  ASUM  E*FFG  J  CONC5300 

C  C3NC5310 

C  CALCULATE    THE    INITIAL    OCwNCOMER    ENTHALPY  C0N05320 

C  C3N05330 

71    HCOC  =  (  (  VFF0+MLLO-Pd3J)*HF+V6eO*H8BO+O.9*(MFF0*ML-L0  )  CONC534U 

6*HV)/(MFF0+MLL0)  C0N05350 

HGHO=HCCO  CONC526G 

C  C0M05370 

C  COMPUTE    THE     INITIAL    3RUM    ENTHALPY  C3N05-330 

C  C0N05290 

FORUMO=HCDU  CON05400 

C  C0N054  10 

C  CONPLTE     NAI  N    BANK    RISER    INLET     ENTHALPY  CON05420 

C  C0N05430 

Mhh0=*LL0  CONC5440 

HJJC=HGHC+Q90/MHH0  CON05450 

C  C0N05460 

C  CALCULATE   TFE    RISER    INLET    DENSITY  CONC5470 

C  •  C0NC5460 

RHCDDO^H^F  CON05490 

~RHCJJO=PHCF  CONC55C0 

C  C0N05510 

C  COMPUTE   THE    C0*NC0ME3    DENSITY  CONC5520 

C  COM C 55 30 

VCCC=(  (MFFC+MLL0-M}BC)*VF+MB3O*V3BO)  C0N055  4O 

£/  (MFFO+Ml  LO  )  CON05550 

RHOCOO=1.Q/VCDO  C0N05560 

RHOGh0  =  5  hOCOO  C0M05570 

C  ■  CONC5  5E0 

C        \  CALCULATE    THE    RISER    OUTLET    CUALITY  CONC5590 

r  C0N05600 
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FILE:    CONSTANT    FORTRAN       PI  NAVAL    POSTGRADUATE   SCHCPL 

HDCO=HCC0  CONC5610 

XFFO^iaiC+fFFO*  (HF-HC03  )  )  f  (**FFC*HFG)  CON056  2U 

XLLO=(Q;?0*MLl  J*(  HF-HJJO)  )  /(MLL  CHFGJ  CON056  30 

C  CONC564C 

C                             CALCULATE    THE    RISER    OUTLET    DENSITY  CONC5650 

C                                                       •  C0NG5660 

RHCFFC»PFCF-XFFO*RhnFG  CON  C  56  7 C 

RH0LLC*RH3F-XLL0*RH0FG  C0N056  80 

C  CONQ5690 

C                             CALCLLATfc    THE    NCN30ILING    LENGTH    OF    THE    RISERS  CONC57GO 

C  CO NO 5  710 

ZBENCi=FNORM-HHD?  CONG57  20 

L0F=ZBENC1  CON0572O 

LDE0=LGF  *A^AXl(lHF-HCDO),0.0)/{  (  HOD  0+  XFF  0*HFG)  -HDOO  )  C0N05  74O 

ZBfcNC2=HNCPM-h^T =CM  CONC575C 

LJL=ZPENC2  C0N0576O 

'     LJKO  =  LJL*\MAXL(  (HF-HJJO)  ,0.0)  /((  HJ  J  O  +  XLL  0*HF  3  )  -H  J  J  3)  C0NC5 7 70 

C                                                                                             -  CONC5760 

C                             CALCULATE    BOILING    LENGTH    CF    RISERS  C0N05790 

C  CON05800 

LEFO=LUF-LUEO  C0N058  10 

LKL0=LJL-LJKQ  C0N35820 

C  C3NC53  2D 

C                             CALCULATE    THE    AVERAGE    DENSITY     INThE    SCREEN    RISER  C3NG5840 

C  C0N058  50 

PhQCFO=(1.0/LOF)*(LEFO/( (XFFO) *VFG)*AL1G(( (XFFC  C0NC5E6C 

6)/\FJ*\(F3+1.0)  +  RHCDC0*LDE0J  CON 058  70 

C  CGNC5830 

C                             CALCULATE   THE    AVERAGE    DENS  IT/     IN    THE    MAn     SANK    RISER      '  CDNC56SO 

C  C0N0590O 

RHCJLO=  (  l.G/LJL  )  *(LKLC/(  (XLLOMVFG)  *ALQG(  (  XLLO  C0NC59  1J 

£*\<FG)  MF  +  l.O)  +  knC  »«0*LJKO  )  C0NG59  20 

C                                "  C0N05930 

C                           CALCULATE   THE    EFFECTIVE    HEIGHT    OF    THE   RISERS  C0NG5940 

C  CCINC5950 

ZDEO=Z0F-LEFO  C0N05960 

ZJM»ZJl-LKL3  CQNC5S70 

ZEF0=ZCF-ZDEG                                                     -  CON  0  59  SO 

ZKL0=ZJL-ZJ<0  CCNG5990 

C  C0NC6CC0 

C                            CCMPUTE    THE    TV.0    PHASE    FLO*    MULTIPLICATION    FACTORS  CON06010 

C                                                                                                                                                .  CON06G20 

PGPAVE*24.794*XFF3**2.0-6.5G66*XFF0+.9776  C0NC6030 

RGRAVK=24.794*XLLC**2.0-e.  5 Jft6* XLLG+ . 9776  CON060hO 

RACLE  =  l5.4564*XFFJ**2.0*lo.4  9'*4*XFFC-.03C07  CO NO 6 C 50 

RACL  J  =  l  5.^564*XLL0**-;'.0+18<.49^-t*XLL0-.00007  CONC0C6O 

RFRICE  =  -J4.0S22*XFF0**2.0«-23.  7ib4*  XFF  0+  -8734  CCN060  70 

RFRICK»-34.08Z2*XLLQ**2.0+23.7164*XLL0+.d734  CONC6C6C 

C  C0N06090 

C                            CALCULATE    SFONO    APPROXIMATION    CF    MASS   FLGw    RATE    AT  CGNC61GJ 

C                            EXIT    3F    SCREEN    R  IS  ERS  M  CON  C6 110 

C  C0NG6120 

PHOEEO=RI-~iODO  C0NQ6130 

*FFCC=<  (PHCC00*G'*Z:C-G*ZDE0*(  (  ^FOODO  ♦  Sh^F  EO  ) /2  .0  )  C0NC614C 

£-G*Z  EFQ*  R  FOE  F_0*R  GRAVE)  /(  (  rCO*  LCD/  DC  J  +  EN"  PCD  C0N36150 

S*eENCCC*EX  nCC)/  U.O  <Ar3*^.J*'<Mr]  G)  +  (  {oh  1EE0-RH0DDC)  CONC6  16  0 

6/CRH0EE0*RHCOOO*A0F**2.0  ) ) +<4 .D*FJE*LOEO *2 .0  C3N0617Q 

-  £  )/  {2.0*mF*{?H~l  =  EC*c  H1DD  0)  *?OF**2.  0)  +CACLE/  CON06180 

£  (RHCEE0**CF**2.J  )  +  (4  .0*F  EF*L  EF0*R  F<i  ICE)/  CON  C6  ISC 

£(  2.0*OOF*RHOEEO*mGF**2.0)  )  )**}.:>  C0NQ6200 

C  C0NC6213 

C                            CALCULATE    SEC~>NC    AP  FR  CX  IMat  [  r>^    CF    MASS    FLOW    RATE  CONCfc223 

C                             AT    EXIT    OF    MAIN    BANK    RISERS  C0N06230 

C  C0N06240 

RH0KK0  =  PH0JJ0  CON06250 

MLL00=(  (  RhjGHO^G^ZGH-^ZJ^O-C  (  RH3 J J O+^HOKKO)  /2.  0)  C0N0626D 

£-G*ZKL0*Ph'JKKJ*Rc3AV!<  )/  I  (  F  Gh*L  O/  OGH  +  EN  T^  GH  COMC62  7C 

x      £*8EN0GH+EXl  1uH)  /(  2.  0*AGn*  *  2.  0*  RH  Q3H0  )  +  (  (  FHGKKO -RHO  J  J3  )  C0N0623J 

£/  (PHCKKJ  »^h"»JJ0»CJL**2.0))+{  4.G*F  JK*LJKG*2.  0  C0NCs2  90 

£  )  /(2.C*DJL*  (RhJKK3*RFQJX)»AJL<*2,0J  -t-^ACLJ/  CQNC63C0 
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FILE:    CCNSTANT    FORTRAN 


NAVAL    POSTGRADUATE    SCHCOL 


C 

c 
c 
c 
c 


c 
c 
c 

c 
c 
c 

c 
c 
c 

c 
c 
c 


c 
c 
c 
c 

c 
c 
c 


6(RHOKKC*AJL**2.  0)  +{  <*.  C*F  KL*  LKLO*RF  R  I  CK)  / 
£  U.O«CJK»Rh")K>i\0*AJL  **2.G  ))  )**0-5 

CDMPARF    PREVIOUS    APPROXIMATION     FCR       RISER     MASS 

FLO*    RATE   TC    CURRENT,     IF    WITHIN    E°90P    CRITERIA    CONTINUE, 

IF    NOT,    UPDATE    AND   ThEN    REITERATE- 

CHECKO.O 

IF  (ABSl  MFFu-MFFCC)  .L1..01  JGO    TC   52 

MFF0=  HFFGJ--'FFO  J/2.G+MFF0 

CHECK=l .C 
62    IFCABS(MLLC-MLL0C).L7..0l)G0    TC   64 

NLLO=  (MLL  )0-lLLG  )/2  .0*MLLO 

GO     10    71 
6*    IF(CFECK.EQ.L.O)GC    TO    71 

COMPUTE     INITIAL    MASS    OF    LIQUID.    IN    DRUM 
DPA5LO*(\.CLD<5M*aHCCCO)/2  .0 

COMPUTE    TNE     INITIAL    DRUM    "LIQUID"    VOLUME 
0M0VC=V0LDRM/2.0 

CCMPLTE    THE     INITIAL    ENERGY   STORED    IN    DRUM    LIOUIC 
CMChlO=C^^SLO*hCPUMO 

EQUATE    INITIAL    FLOW    RATES 

MCDO=MFFC 
MGhO="LLQ 


6C0 


6Ci 


CQMPLTE 
CRUM 


INITIAL    .MASS    CF    STEAM    IN    STEAM 


0STM0=VDLQRM*RHOV/2.  C 

CQMPLTE    HXFFO    AND    HXLLO 

hXFFO=HF*hFG*XFFO/Z.O 
HXLLC-HF+HFG-XLLC/2.  C 


WRITE 
FORMA 
WRITE 
WRITE 
FHPM  \ 
WRIT  E 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRI  TE 
WR  ITE 
WRITE 
WRITE 
WRITE 
WRITE 
WR  ITE 
WRITE 
WRITE 
WR  IT  E 
WRITE 
WRITE 
WRITE 


IS  ,60 
7(  1HI 
(o,  IN 
(6,60 
T(  IHQ 
(6,  M 
(6,60 
(o,  I  1 
(6,60 
(6, CO 
(o,  60 
(6,CC 
(6,  6C 
(6, CI 
(6  ,60 
(o,C0 
(o  ,60 
(  o  fC  C 
(0,60 
(6,60 
(  fct'lU 
(7,  IN 
(  7, IN 
(  7,  IN 
(7,CC 
(  7  »CC 
(7,  CO 


0) 

) 

CON  1) 

L) 

) 

CT'l^J 
I  J 

CON  3) 
1) 

4STL) 
1) 

NST2) 
I  ) 

N5""3  i 
1) 

NST4J 
I  ) 
N5T5) 

1) 
0) 

"PUT) 
CON  1  ) 
CCN21 
C-3*l  I) 

NST1J 
NST2) 
NSTji 


C0N063  10 
CON06320 
CON  06 3 30 
C0N06340 
C0NC6350 
CDN06360 
C0N063  70 
CJNC63E0 
COM  06  390 
C0NC64  00 
CDNC641G 
CO  NOo 4 20 
CON06430 
CON0644Q 
C0N06450 
CONC6460 
COM06470 
C0NC6480 
CONC64S0 
C0N0t>500 
CON0651D 
C0N06520 
CGN06530 
C0NG6540 
CON  065 50 
C0N06560 
C3NC6570 
C0NC6580 
C0N06590 
C0NC66CC 
C0N06610 
CON06620 
COMC6620 
C0N066A0 
C0M06650 
C0N06660 
CO NO 6 6 70 
C0N066E0 
C0NG6690 
C0N0&7  00 
CONC67L0 
CON06720 
C0NCo730 
CONC6740 
C0NC6750 
CON06760 
CON C6 7 70 
CON 0  67 80 
C0NC6790 
C0NC68C0 
C0N068  10 
C0NC6E2C 
C0MC68  30 
C0NCt>840 
C0NC6650 
C0N06860 
C0NG68  70 
C0NGO880 
C0N06  890 
CON C6 SCO 
C0NC6910 
C0N06920 
C0NC6S3G 
CQN069A0 
C0NG6950 
C0NC6S60 
C0N06970 
C0NC6SEC 
C0NC6990 
C0N37000 
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FILE:    CONSTANT    FORTRAN      PI 


NAVAL    POST5R4CLATE    SOGCL 


WRITE(7,CCN'T4J 

WRITE  (  7.CCNST5) 

WPITE^.^UT^UT) 

STCP 

END 


CONC7C1C 
CQNC7020 
CQN0703O 

CDNC70AU 
CON0705J 


APPENDIX 


B 


C6/2C/1*       21. 26. IS 


FILE:    CS  f-F 


FORTRAN       PI 


NAVAL    »0STGR4JLATE    SCHOOL 


//HOPEFULL    JH    I  2473,C5e7,\E  31  ,,9)  ,•  WALKER    ^'<c    1319',TI^E=2 
//       EXEC    C5^'°XV 
//*.  SY5IN      DO       * 

* 
* 
* 

* 

*  THIS    0    TY°E    81ILER    M33EL     IS    W3  I  TTCN    IM   CS^P-III 

*  CC*FJTE«    LANG'JACfc..      IT     IS    A    \j-»m  l  I'JE  ?"    M1DEL 

*  WHICH    UTILIZES    A     SfcT    jf     INITIAL    CC\'OIxI  QNS     ANC 

*  CTHEP     3AT4    PREP  fi5£D    3Y    VI     Initial    C  CNOI  TI  ON    FORTRAN 

*  NCCLL    ViMCh    PRECEDE*    THIS    ONE.       TdE    M  JOEL     IS     GQOO 

*  FOR     S^iLL     TFA^SIENYS     AF-JsjNO    TrtE     CPfcRATiNG     POINT 

*  SPECIFIED    BY    THE     INITIAL    CONDITION    PROGRAM. 
* 

INITIAL 

rONST    DR  A4l  =  C.  C  ,DRA42  =  L.  3,uFA4  3  =  0.  3  ,  CRAt4  =  0.  )  , ORA45=0.0,  ER*46=0.0,  ... 
OR  /4 7=0.0,  C? £4d=u  .0,  0^  £2i=C.0»DRA32  =  C.  Ct  OR  A  33-0.  0,0*  A  34=  C.  C,... 
§RA3  5  =  C.  C  ,CPA36=J.0,CR/37*0.0  ,CR  A38=0.0 
CONST   0NE=1.  C,     ... 
TRRO=      2563.0000  ,TV7J=      oOO.OjOCO  ,  TWWO*       1033. OOOC  ,TYYC=.. 

515.4274S  ,TZZO=      600.99976  ,TXXO=       360.3o450 

CCNST    TW0=2.0,     ... 

DA<JHO  =       ?"i54667.0  ,CMASSO  =       5415.3164  ,,jiaAJ=       45.00C000    ,TAAO=.. 

254.CCC0b  ,*FQQO=       3.307't999  fMAOOO=       55.398380  ,MNNO=... 

7.7500000  rMMMO=        ^tb. 000000  ,  V*LVrO=     0. 505V9959  ,MQQO=... 

55.2C6375  ,CVN  =    0.&5523023  ,CA3=       1  .0  76  WOO  ,  C*V  =  0  .  1099999  5  ,..  . 

C'Ys    0.10999995  ,r.?x=     C.1CS9S995  ,CPZ*    0.10999995  ,C°^  = 

0.  10999995  t  D*0HL0=396i62o.0i  D'*  AS LO  =300 2  .2000, DMOVO=  163.65698 
CONST    TH°EE=2.C,     •  •• 

DROCFO=    0.0  ,  JcnJLO=    0.0  ,MM-1E=    0.  449S  <=979E-G2  ,  .  .  . 

*Gh3  =      662.70825  ,MR=>3  =      ^'i.OJOJOJ  tlSnry^ 

lc.c6<SS£  »OSlMDMa       60.C00300  , «CD0=980. 02 734, CLCV0=163 .66 

CCNST    FH'ift  =  4  .0,      ... 

KSJL=     C.5<99SS91E-04,l<.SGF  =    3  .2^99  9996  E-04  ,  KS  Gh=    0  .49999999  E-05  »  KSCD=  .. 
G.49S9S999E-05,KSJE-=     C.  2  c.  S99996L  -0  A  ,  nSE  F=    0.  299S99  96 1-04  ,       ... 
PSJK=0.59999791E-D4,     ... 
KSKL=     C.5<SSS«S1E-04,AJI.=       10.678426  ,0JL=       19^.57166  ,LJL2=... 

23420.000  ,LJL1=       IG.CGOOjO  ,30F  =       35.  897476,  LOF-15.  75000C  ,..  . 

LDF  1=       12.7CGC30  ,iCF-      t. 0646043  , LGH=       25.00JOOO  ,AGH=... 

1.4451332,SAMOL,=  250,RIS7IM«10.0 
CCNST    HIVF=5 .0  ,     ... 

OGH=        1.9183321  ,0C0=       4.6t>99991  ,LCC=       9-*. 503330  ,LDEO= 

2.6166D15  ,LJI\0=     C.6C5C0102  ,ACD  =       2  .  95  47  77  3  ,  MA  S  S  V=  7951  .  3  984 


,VASSW: 


lei  i.  CCCC 


13752.000 


, VASSY= 


^192  .000 


S  jX  = 


,Z3E^C1=       15.750C30 


,f-C3=    0.8  1J23693E-01  , 


24442. GOO, '4ASSZ=.  . 


FCE  =    0  .92714846E-01,  FEF=     G  .927  1  4a4tjE-  01  ,F  GH  =  G.  8  1023  t  93E-C  1 
CCNST    SIX =6.0,     ... 
ZBENC2=        12.41CCC1  ,HA3SH=    71.3S5  ,VOLJL=       1 06. 79426  ,VC  LC . 

F=    51.6234-68  rvrt_DRv  =        227.31396  ,VHLHJ=        93.174744 

CONST    £EVEN  =  7.0,    ... 
CQ3  =    J. 31595814  ,CFS=     C.2045367C         '    ,Z  ST=     C. 29342304  ,CTU  = 

C.2602606^,  ,CNP=    0.5  3347189  ,KMP=     13.508d20  ,  K  W=  2.  099772  4  ,  .  , 

SIGMAA=     C.  21122S12E-CS,<sRS  =    G. 87139436  ,KST  =      2.2037563 

CCNST    fciG»-T=3.0,     ... 

KTU=    C.441266CC  ,£NTPGF=   0 .99999964E-0 1, ENTRCD=     0. 9999 SS54E- 0  1 ,  .  .  . 

8tN3:0=    C.S9999964fc-01,EXITGH=       1 . 0 0000  CO , E < I TC 0= ..  . 
1.0030300  ,  HA40=        221.0J0JJ  ,KTV1=     0.  1064^3  7oF- C  1  , .  .  . 

KJN4=    .1C777515  ,F  JK>=    0.981  96  149E-31  ,  FKL=    0  .93  Wo  l-»9  E-01 ,  .  .  . 

GC=        32. 199V97  ,TAM8=        80.  CCCOGu  ,ZCD=       15.753000  ,ZDF  = 

15.753300  ,ZJL=        12.410001  ,LJL=        12.410001  ,... 

^ENJGH=C.  149i9SSd,NCN3  =  t'«7.714  60,G=3  2.l999  97,0STv3=2  43.59  8O2 
CCNST     NP!E=9.0,     ... 

KX=        IG.  077373  ,KZ=  14.930669  ,  KV=    J  .  32  769o  1  3F- 0  5  ,..  . 

FHV=       192CC.4S6  ,ODE=       35.89747o  ,OKL=       196.57166  ,  DEF*  •  • . 

35.897476  ,  DJK=        19o  .5  7  1  b6  ,i-»  Y=  0.  26682  50  2E- 02,  .  .  . 

ORA1C=C.C,ORA2  0  =  0.0,XJ:F9  =3  .13  9503  30  E-01  ,  ... 


CS'<3O010 
CSVC0020 
CSMCCC20 
CSM00040 
CSWCCC5G 
CS'00060 
CS  '-'000  70 
CS^CCCSO 
CSMoC090 
CSN'00100 
C  SMC  CI  10 
CSV"00120 
CSV00130 
CSM00140 
CS^00150 
CSMCC16C 
CS"00170 
CSN00180 
CSMCC190 
CSviCC2C0 
CSV00210 
CSMCC220 
CSM30230 

•CSV30240 
CSMC0250 
CS*C0260 

.C  SMCC27C 
CS  ^00280 
CS^CC290 
CSVCG3C0 

.CSV00310 
CSMCC32C 
CSM00330 
CSVC0340 
CSMCC3  50 
CSV00360 
CS^00370 

.CSM0038O 
CS  ^00390 
CSVCG400 
CS^004  10 
CSV00420 
C  SVC0430 
CS  ^30440 
CSVC0450 

.CS'iGOAtO 

.  CSK004  70- 

.CSMCC4E0 

csyccsco 

CS^COSIO 
.CS^30520 
CSWCC530 
CSN'00540 
.CSN:C0550 
•CSMCC560 
CS. 400570 
CSV00580 
CSMCC5SC 
CS^00600 
CSNC06  10 
CSMC0620 
.CS^00630 
CSMCC64C 
CS^C0650 
CSWCG660 
C5MCC670 
CSV00680 
C3^C0690 
CSM00700 
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F  I  LB'.    CS*P 


FORTRAN      PI 


NAVAL     POSTGRADUATE     SCHOOL 


CCMPUTE   Th£   TOTAL    FL'JE    GAS    FL^W    RATE    I^TD    BOILER 


COMPUTE    THE    ENERGY    ENTERING    THE    BOILER 


CONFUTE   THE    ENERGY    TRANSFERRED    TO    THE    SCREEN    RISERS 
VIA    RAOIATICN 


XL  L  G=  C.4CCCCCCCE-Cl,MFFC=SiC.  02  7  34,  VIL0=662.  70  825 

CCf-ST    TEN=10 .3  .... 

"5RAll=C.  C»DR A  12=0.0  ,ORA13  =  C.C,DRA14=C.O»DP/15=O.0»CRAi6=O.0t«... 

DsA  17  =  0.0  ,0"  *2l=  C.  u,^:  A2  2=C.0,")R  A  23  =  0.  0,0K  A2  4  =  0.  0  t.  .  . 

FCL=«015»CC0=63836.3  95,ZGF=12.4100Ol  ,  ... 

0^A2f=C.  C»0=  A2fc=C.  0»  OR  A2f=  CO, H<LLC  =  505.  31  ,HXF  F0  =  502.  34 

CYNAMK 

* 

*  INPUT    EQUATIONS 
* 

VALVE=.  015*R£VP(1D.C)-.015*RAMP(20.0  J+.51 

«<\»N=MNNC 

v/A=MAAO 

1AA=TAA0 

* 
* 
KC=MCCO 

* 

4 

* 
00=000 

* 
* 
C1=SIGMA*  UTRP+460.0  )**4.0-{TVV  +460  .0  1**4.0) 

*  COMPOTE    THE    "ATE    EQUATION    FOR    FLRNACE    FLUE    GAS 

*  TEMPERATURE 

* 

DT  PR=  (CQ-C1-'1R«*C0R*(TPR-TAM3  J  )/(MASS0R*20R) 

»  CT-1PITE    THE    FLUE    GAS    TEMPERATURE 

* 

TRRMNTCRLI  7RR0.DTRR) 

4> 

*  CCMF'.TE    Tl-£    rE^FERATJRE    OF   THE    FLUE    GAS 

*  LEAVING     THE     SUPERHEATER 
4 

PHIi=?.C*VPR**C.4*C«S/KRS 
TSS=  (TRC*(PHI  1-  1.  0)  +2.  C*TU**  I  /(P-Hl  +  1.0) 
* 

4 
* 

Q3=MRR*CRS*(  1RR-TSS) 
* 

*  CCVFUTE    Thb     » A I  N    BANK    EN  ERG*    TRANSFER 

*  FLUE    GAS    TJ     TUBE    METAL 

* 

Q5  =  MSS*C  !!*<  TSS-TTT) 
* 

*  CO^o-jTE    THE    TEMPERATURE    OF    THE    FLUE    GAS 

*  LEAVING    THt     VAIN    BANK 
4 

CHI2  =2.0*'1SS*40  ,4*CST/KST 

T7T=(  TSS-M  PHI2-1.0)  +  2.0*TYY)/  (PHI2  +  1  .0) 

*  CCVFUTE    THE    TEMFEPATURE    OF   THE    FL'JE    GAS    LEAVING- 

*  THE    ECONOMIZE* 

* 

PH13=2.C*C  U/K1L 

TUU=  (TTT*(°HI3-1.C)*2.  C*TX>)  /(P4I3+  1.0) 

*  COVPLTE    ECCNC^I  2ER    ENERGY    TRANSFER 
4  FLUE    GAS    T)    T'JPE    METAL 
* 

Q7=MTT*C  1LMTTT-TUU) 
*) 


CCyFL'TE    ThE     SUPERhEATER    ENERGY    TRANSFER 
FLUE    GAS    TO    TUBE    META  L 


CS^00710 
CS^CC720 
CSM0O730 
CS  *CC740 
CSMCC750 
CSM00760 
CS"JG770 
CSM00  780 
CSM00790 
CSVC0300 
CSM00810 
CSN.003  20 
CS"CC630 
CSVQ0340 
CSVC0850 
"  SMCC660 
CS '4  008  70 
CS"CG880 
C  SV 00890 
CSV00900 
CS^OCSIO 
CSMCC920 
CSM009  30 
CSVCCS40 
CSM0O950 
C5  VC0960 
CSMCCS70 
CSM00980 
CS"0C990 
CSMC1000 
CSVGIOIO 
CSMC1020 
CSMC1C30 
CS"01040 
C  SVC1C50 
C3«0  10  60 
CS  f-:iG70 

csmciceo 

CSVC1090 
CS NO  11 00 
C  S^  CI  110 
C5"'Jl  120 
CS*C1130 
CSMC  1140 
CS''G115J 

cs*cn  to 

CS^'0  1170 
CS "CI  180 
CSMCl  l?0 
CS'-'012O0 
C  S>-G1210 
CS.vO  12  20 
C* "C1230 
CS^C1240 
C SMC  12 tO 
CSVJ1260 
CS^C1270 
CS''G12  8Q 
CS  N CI  2 90 
CS»-C13C0 
CS^OIJIO 
C  <VG1320 
C  SM01330 
CS*'01340 
CSMC1350 
CS'-1013CO 
CS«C1370 

c$mci3ec 

CSM  0  1390 
CSVG1400 
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FILE:    CS"d 

* 
* 
* 


FORTRAN       PI 


NAVAL    POSTGRACUATE   SChOOl 


ISE    CONT 
FLOW    =AT 


*L 
("f 
OR 
LP 
PS 
IF 
09 

DR 

oa 

OR 

C 

CR 

DR 

CR 

OR 

OR 

CR 

DR 

CR 

OR 

DR 

CR 

Ok 

OR 

OR 

1 

fcN 

FR 

IF 

CR 

CR 

CR 

DR 

CR 

Ok 

DR 

CR 

OK 

CR 

CK 

DR 

Ce 

DR 

DR 

DR 

DR 

Dp 

CR 

2 

EN 

* 

* 

RM 
RH 


L=MJJ- 

f=*CC- 
Al 6=0£ 
A2(:=CE 
CCECIR 
(KEEP. 
FDJL=( 

ORA 
4  50=DR 
£52=CR 
A51=CR 
A3  1=CR 
£3  2=  CR 
A33=DR 
A34=CR 
£3  5=0R 
A3t=DR 
£J7=CR 
A3£=DR 
£11=  CR 
£12  =  CR 
A13=DR 
£14= CP 
A15  =  0R 
Alfc=0R 
Al  7=CR 
CCNT 
C'RHCE 
CCfcCUR 
(KEEP. 
F3CF=( 

ORA 
663=C  = 
£6  2=  CR 
A61=DR 
£4l=C? 
£42  =  CR 
A43=0R 
£4<«=  CR 
A45=0R 
£4  6=  C 
£4  7=CR 
A4E=DR 
A2  1=  CR 
A22=0R 
A23=DR 
A2  4=CR 
A25=UK 
A2*=DR 
£2  7=CR 

COW 
CPRQCt 


DRHOJL 
C^FnCF 

pi  wj; 

?  1V(  DR 

t     CPFC 

NE.l  )„ 

D?  M-J  + 

3H-CRA 

A52 

A51 

£31 

^32 

£33 

A34 

\35 

£36 

A37 

£38 

All 

*12 

A13 

A14 

A15 

A16 

A17 

A18 

INIE 

DU»E 

£     CPFC 

NE.  1)  i 

CR^S  + 

4l+C*£ 

A62 

A61 

A41 

A42 

£43 

A44 

£45 

A4  6 

A47 

£4  8 

A21 

A22 

A2  3 

A24 

A2  5 

£26 

A27 

£2  8 

TINLE 

CURE 


*VOLJ 
*VnLC 
£10, R 
A20,? 
JL  =  M 
C  TO 
CR  \  17 
J2+C3 


LJUITY     RELATICNSH1F   TG    COMPUTE   THE 
E    UP    THE    5  ISE&S 

L 
F 

H3JLJ 

HOOF  ) 

LTP1  (DRU4) 

1 

+  DR  I-  16  +  0=  M5  +  DP-  A14  +  D*  A13  +  3P  A12+3RA1 

A33-frC3A34+CP  A35+CRA3t»+0RA37+D=A  3d  )/ 


MASS 


16.0 


CF=  F  I 

a  td 

CkA27 

42  +  C3 


L'R^(DR£24) 
2 

+  3?  u^o  +  C°  A25  +  D:i  A24+D3A  2 3+D  R A  12+  DR  A2 
£43  +  CR  £4^  +DRA45  +D^A4y  +DP£47+DRA43 )/ 


1+.  .  . 
16.0 


COMPUTE    THE     £VEC£GE    D£NSlTY    IN    THE    RISERS 


3JL=(  1.0/LJL>*( (L 

+ShCJJ*L JK) 
C3F=(1  .9/LCF)*  (  (L 

+^H1CC*L0E J 
CFF==FCf-XFF«PhlF 
CLL*RHOF-XtL*RHOF 
COMPUTE 
TUBE    MET 


KL/IXLL*VFG)  )*  ALOG  <(  XLL/VF  )*VFG  +  l  .0 
EF/(XFF*VFG))*ALOGUXFF/VF  )*VFG*l.O 


I... 


RH 

RH 

* 
* 

C6  =  KY*PS£T**U.O/2.C)*(TirY-TSAT)  **3.G 


G 

G 

THF    MAIM    RANK    ENc^GY    TRANSFER    - 

AL    T  C    W£IN    BANK    M  1XTU"  E 


CSK141C 
CS"  014  20 
CSVC1430 
CSMC1440 
CS MO  1450 
C SWC1460 
C3V'014  70 
CS  V0  1480 
CSMC14S0 
CSMC15G0 
CSVJ1510 
CSMC1520 
CS  "0  1530 
CSM01540 
: SMC  15  50 
CSVG1560 
CSV01570 
CS**01530 
CSVC1590 

cs'<cieco 

CSM016  10 
CSW01620 
CSMC1630 
CS  "01640 
CSVC16  50 
Cb,4Cl660 
CSN'ul6  70 
:SMC1680 
C.S"J16  90 
CS"C1700 
r.  SMC  17  10 
CS  wo  17  20 
CSV  CI  730 
C5MG1740 
CSVC1750 
: S^C  I  760 
CS NO  17  70 
CS"C1  780 
C3"01790 
CS  NJ1800 

cs^cieio 

CS1'' -01820 
CSVG18  30 
C  S^CIBAO 
CS  "01850 
CSvO  i860 
C S M !  G  1 8  70 
CS  »» j  ij  80 
CS"CIPS0 
CV019  00 
CS"U19  10 
C  S°.Clf»20 
CSV01930 
CS^OIS^O 
;SU0  19  50 
CS  VJ19f>0 
CS''CIS70 

cs^oiseo 

C^ V0  19  90 
CS-'C2CC0 
CS" 020  10 
CSVC2  020 
CSMC2C?0 
CSV  Vjz,j 
C  SH020SO 
C%v020t0 
CSv02  0/0 
CSM-L2060 
CSMC2090 

csvo^ioo 
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FILE:    CSMP 


FORTRAN      »l 


\AVAL    PC3TGR4CL/TE    SCHOOL 


CCM0UTE    THE     SCREEN    RISE*    ENERGY    TRANSFER    - 
TUEE     1ETAL    TO    SCREEN    RISER    MIXTURE 


* 
* 

02=KV*PSAT**{  4.C/3.CJ*(  7VV-7SAT)  **  3.  C 

* 

*  CJ*PITE    THE     FISER    AVERAGE    ENTHALPYS 
* 

hDF=H0D*XFF*HFG/2.0 
I- JL=HJJ+XLL*HFG/2.  0 

* 

*  CCMPLTE    THE    RISER    OUTLET     ENTHALPYS 

* 

hFF=FF+X  FF*FFG 
hLl=HF*XLL*HFG 
# 

*  CONFUTE   THE    RISER    EFFECTIVE    HEIGHTS 
* 

<CE=ZCF-LEF 
ZjK=ZJI-LKL 
ZKL=ZJL-ZJK 
ZEF=ZCF-ZCE 

*  COMPUTE    1HE     SPECIFIC     VOLUME    OF     THE    ECONOMIZER 

*  CJTLE^    LIQUID 

* 

VRB=.0160C4oa-.CGC0G20I<i6*T3B+.O0O00C036  511*"PB**2..O... 
-8.i42E-ll*Tee**i.)  +  i..i0iE-lJ=s'T3B**^  .0- 1  .  143E- lb*T3B  **5.  0.  .  . 

*  c.034E-2C*T33**o.  G 

*  COFITE    THE    RAXE    fcOIUTICN    FOR    RISEQ    OUTLET    (DUALITY 

* 

ChXLL  =  (MGG*(hJJ-FF-XLL  *hFG/2.0  J+Q6-  MLL*XLL  *HFG  /2.  0)  / (RHO JL* VOLKL ) 

riXLL  =  I  NTGRLIHXLLCtChXLLJ 

aLL=( H>LL-FF |*i.U/HFG 

CHXFF=(MCC*(hCC-FF-XFF*FFG/2.0  )♦  0<i-M  FF*X  FF*  HFG/2  .0  )  /  (  R  HOOF  *  VDLEF  ) 

HxFF=!  NTGRH  HXFFG.DHXFFJ 

XFF=  (FXFF-FF J*2.D/HFG 


CCMPLTE  ThE  RISER  BOILING  VOLUME 


* 
* 

» 

VCLEF=VCLCF*LEF/LCF 

V3LKL=VJLwL*;tKL/LJL 

A 

*  COvF:!TE    THE    NCNEPILING   LENGTH   OF    THE    RISERS 

tDE=LCF*IhF-FCCJ/( (  HF  +  a  F  F  *  FFG  )- HDD ) 
lwK=LJL*  (HF-rlJJ)  /((hF  +  ALL*»-FG)-HJJ) 
* 

*  COM°yTE    ThE    RULING   LENGTH   3F    THE    RISERS 

d 

LKL=LJL-LJK 

IEF=LCF-LCE 

* 

4 

* 

*  COMPUTE    THE    MASS   "ATE    EQUATION    FOR     STEAM   CONDENSING    IN 

*  ThE    CRUM 
* 

f»CCND=560.93*(PMV/(TMM*460  .0  )  **J  .5-°  SAT/  (  TSiT+46D.O)**G.5J... 

+.C2568 
* 

*  CCMF'IT£    THE    RATE    EQUATICN    FOR    DRUM    L  IQ'J  ID    MiSS 

COMAS  L=MLL*<  l.J-XLL  }-MFF*{  1  .C-XFF  )  +MCOND+-M  B3  0-MCC-  MGG 

*  COMPUTE    THE    DRU*    LIOUIO    MASS 


CSMC2110 
CS^02120 
CS*02i30 
C  SM02  140 
CS  "02150 
CSM021C0 
CSMC2170 
CSM02180 
CSMC21S0 
CSM0220Q 
05^02210 
CSMC2220 
CS* 022 30 
CS*0<!240 
CSMG2250 
CS  ^02260 
CSVC2270 
CSMC2280 
CS^C2290 
CSMC2  3C0 
CSM02310 
CSMC2320 
CSMC2330 
CS'<02340 
CSVQ2350 
CS*G2360 
CS  "023  70 
CSMC23  60 
CSVC2350 
CSVC2400 
:S''C241G 
CSN' 02420 
CSNC2430 
CSMC2A<.0 
CSW02450 
CS "C246C 
CSM02A70 
CS^C2430 
C  SMC24S0 
CS "02500 
CS"J2510 
C  SMC 2 520 
C5 "0  2  5  30 
CSVC2540 
CSMG2550 
CS "02560 
C  SMC25  7C 
CS,-'02  5eO 
CSVC2590 
CS'-' C26G0 
CS"Q26  10 
C  Sv0^620 
CS •■'  J2630 
CS "D2  640 
CS"C26  50 
CS'"C2t60 
CS '026 70 
C  S"  C2fcEC 
CS"C«i690 
CS  "ZZ1CQ 
CS--C2710 
CS" J2  7  20 
C 3 ^02 730 
C 3" 0^740 
CS  "02750 
C  S"C27fcO 
CS"C2770 
CS''027B0 
C  SMC27S0 
CSMU2800 
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FILE*    CSMP  FCRTRAN       PI  NAVAL    PDSTGRACLATE    SCHOOL 

j*.ASL=lNTG3L(0yASLU  ,00 MAS  L) 

*  COMPUTE    TH£    PATE    EQUATION    F0=>     ENERGY    IN    THE    D-^U*    LIQUID 

CJuOHL="LL*(  l.J-XLL  J  *HF*  IFF  *(  I.  O-XFF  )  *HF+MC0N0*HFG+M3B*Hq3.. 
-KC*HCC-»'%iij*HGG 

*  COMPUTE    THE    ENERGY     IN    THE    DRUM    LIQUID 


CMCHL=  INTG»L (DMOHLOfDOMDHL) 


COMPUTE    THE    ENTHALPY    CF    THE    DRU*    LIQUIO 


* 
DHsOPOHl/CJ'JSL 

* 
* 

*  CJMPITE  THE  DRUf  SPECIFIC  VOLUME 
* 

CRYUND=PCU* <M.L*VFF) 

PPOCEDLRE    RISE2=FI  LTR4ICRYUND) 

RlSEi=0EL4Y{  <:50,RISTIM,CRYIjND) 

RI  S£2=CRrUf>0 

IF(  VALVE.  G7..  5DG0    TO    55 

PlSc  =  RIS  E2 

GC    TC    57 

55     -  RISE  =  RISE1 

57  CCMINLE 

ENC°9nCEDJ^  E 

CCOV=  (  <  'Ff-*VLL-  184)*VF*5CI»(  vFF+MLL  )  *VV  *MBB  0*VPB-  (  MCC*  MGG)  *  VF  *MC  UNO  . 

*  VPG-»I  SE*VV) 
CM0V=IN'TGRL(O  1D70,D0M0V) 
* 
* 

LEVEL=(0*»DV-VCLls.v/2.0  M  t  LST  ^0M*CSTMDM  ) 

*  COMPUTE    THt:    03WNCJMER    ENTHALPY 


CJMPLTE    THE    DRU  *    LEVEL 


HCD=I  (  ►FF+PLL-MB8)*HF*ME3*^BB)/  (  NFF  +  ^LL) 

FGH=  FCD 

HCC=HCC 

HCC»HCO 

FGC-HGH 


* 

V*Q=t  (  NFF  +  NLL-I"BB)  -VF+NE3*  V86)  /  (  ^FF+WLL) 

°hizc=  l.o/vr.o 

RhCGh=RHCCO 

* 


COMPUTE    THE    DQhNCCMER    SPECIFIC    VOLUME    ANO    DENSITY 


COMPUTE    THE     SATU^.ATIDN    P^ESSU'E     \  NO    TEMPERATURE 
CCPFfcS  PONCING   TC   ThE    DOWNCOMER     ENTHALPY 


PS  AT=£XP  U*L0C-<FSAT)-4  .hj7C3  )/.  26452) 
TSAT  =  E  >Pt  I. 221 51* ALOJI PS AT) +4. 77 123  J  ) 


* 
* 

4 

HJJ=HGH*C9/MH 

* 

4 


CCMF'ITE  Tht    EJThALPY    CF    ThE   LIQUIO    ENTERING    THE 

MAI N    BANK    RI SER 


CCFUTE   ThE    RATE    EQUATION    FDR    THE   MAIN    BANK    AND 
SCREEN    RISER     fUEE    METAL    TEMPERATURES 


DTVV  =  (C1-C?)/  (  fASSV*CFVJ 
0TYY=(Q5-QO  /(  MASSY*CPY) 


CCVPLTE    ThE    SCREEN    ANO    «A1N    BANK    RISER    TUBE 
METAL     TFMPEFAFUKES 


CSMC2810 
CSV02320 
CS V02330 
CSMC2840 
CSMC2650 
CSM02860 
CS*C287G 
CSM02  3  80 
CSMC2990 
CSMC29C0 
CS  "02910 
C  SMC2S20 
C3MC2920 
CSMC2  940 
CS"C2S50 
CS"-'029  60 
CSMC2970 
CSM02980 
CS "02990 
CSVC2CC0 
CSMC3U10 
CS MC3020 
CSMCiC30 
CSMQ3040 
CSM03050 
CSMC2C6C 
CS"C3070 
CSMC3C80 
CS  "33090 
C3*C310O 
CS'-'C31  10 
CS^03120 
C  S^C3l 30 
CSMQ31W 
CS "C3150 
C  SMC  J  160 
CSV03170 
CS"C31S0 
C  S*03193 
CS "03200 
CSMC3210 
C3MC3220 
CS MO  32 30 
CS^C32hO 
CS*'0  22  50 
CSf03260 
Co''C32  70 
CSV  J  3230 
C  S'^C  32  SO 
CSM 0  3300 
CS "Cj3  10 
C SMC  J  320 
Ci"  Jj>330 
CSvCj3^0 
3  5MQ3350 
CS  v^3360 
C  SMC32  70 
CS'' 03380 
CSV Jii90 
:  SMC34C0 
CS"G3A10 
CSMC2420 
CSMC34  20 
CS"C3440 
C  S^C3A50 
CSM0iA60 
CS"C3A70 
CS"o2<.6G 
CSM0  JA90 
CS"C35CC 
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PILE:    CSmP  FORTRAN      ?l  NAVAL    POSTGRADUATE    SCHOOL 

TVV=INTGRL<TVVC»QTW)  CSM03510 

"YY=1NTGRL  (TYYJt  CTYY)  CSM03520 

*  CSMCJ5  20 

*  COMPUTE    THE    HEAT    TRANSFER    FOR    THE    Q  ES'JPERriE  ATER  CS"Oj540 

*  -STEAM    T3    TUBE    NETAL  CSMG3550 

*  CSMC3560 
Q9=CNP*MNN*(TNN-TPP)                                                                                                       .  CSM03570 

*  CSMC3560 

*  CONFUTE    THE    TEMPERATURE    OF    THE    STEA*    LEAVING   THE  CSMC3590 

*  CtSUPEPhEATER  CSWC3600 

*  CSMC3610 
TPP=(TNN-TZZ )/(ExP(KNP / ( CNP *MNN**0.2 ) J ) * TZ2  CS*036  20 

*  CSKO3430 

*  SE1    THE    DOWNCOKER    E^RANCE    AND    EXIT    T  EMP  ERATUR  ES  C5M036A0 

*  EQUAL    TO    THE    SATURATION    TEMPERATURE    CORRESPONDING  CS*03650 

*  TC    CRUf    ENTFJLPY  CSMC36tO 

*  .  CSMC3670 
TGG=TS4T  CSM03680 
1HH=TSAT  CSMC36S0 
TCC=TSAT  CSM03  700 
TDC=TSAT  CSM03710 

*  CSM03720 

*  COMPUTE    THE    E\lEr!GY    TRANSFER   FOP    THE    DESUPERHEJT  ER  CS^J3730 

*  TUEE   METAL    TC    DRUM    LIQUID  CSMC3740 

*  CSM0375O 
Q99=KZ*(TZZ-THh)  CSMJ3760 
»  CSMC3770 

*  COVPUTE    THE     PATE    EQUATION    FOR    DESU  FEP.Ht  AT  ER   TUEE    METAL  CSMG376G 

*  TEMPERATURE  CSV03790 

*  CSMC3  6C0 
CT  ZZ=  (09-Q  99  )/(  MA  SSZ*CPZ)  CS  VQ3810 

*  CSMC2820 

*  COMPUTE  THE  DcSLPERHEATER  TU3E  "E^AL  TEMPERATURE  CSMC3830 
<  CSM33840 
TZZ=INTGRLITZZO,CTZZ)                                           .  CSMC3E50 

*  CSUQ3860 

*  CT^PUTE    ThE    ECONOMIZER    ENERGY    TRANSFER-TUBE    METAL  CSWC3870 

*  TO    FEcu    WATER  CSMC2880 

*  CS'1U3890 
C8  =  M/A*C£E*(TEE-7AA)  '  CS^0i90O 

*  •  CSM039  10 

*  COMPLTF    THE    FEED    TEMPERATURE    AT    CUTLET  CS^J3920 

*  CF     ECONOMIZES  CSVC3930 

*  CS^'C3SA0 
TB8=(TAA-TXX)/<EXP(KX/  (  C  A° i*MA  t  **  C.  2)  ))  +  TXX  CSV039  50 

*  CSM03960 

*  COMPUTE    THE    RATE    EQUATICN    FOR    THE    ECONOMIZER    TLEE  CSv0i9  70 

*  V£T£L    TEMPEP^TUPt  CSMC3980 

*  CSMC3S90 
DT  XX  =  (  08  -Q  7  )  /  (  M  A  S  S  X  *C  P  X  )  CS  *>. 040  00 

*  C SMC AC 10 

*  COMPUTE    THE     ECONOMIZER    TUEE    METAL   TEMPERATURE  CS"J4020 

*  CSM04030 
T>>=INT^PL(TXXC,CTXX)  CSMC4040 

*  CSM0^0  50 

*  CCNF'J-E   THE    SUPERhEATED    STEAM    OUTLET    TEMPERATURE  CSVC^CfcO 

*  •  CS'-'040  70 
tmn=  (TMM-TWW  J/(  EXP(  Kvi/  (CMN  *MNM  **  C.  2)  )  )  +TWW  CS  *C-tC8'J 

*  CSMCA090 

*  CCMPUTE    THE     SUPERHEATER    ENERGY    TRANSFER-TUBE    M£T4L  CS'-'QALUO 

*  TO    STEAM  CS"C4llO 

ca=c  *n**m»*  hnn-tmm)  c  smos12o 

*  CS"04130 

*  COMPTE    THE    RATE    EQUATION    FOR    SUPERHEATED     TUBE    METAL  CS*C*140 

*  TEMPERATURE  CSMC4150 
4  CS*-'04l60 
CTWW=(C3-Q^)/( NASSW*CPW)  CSWC417C 

*  cs^o-^ieo 

*  COMPUTE    THE     SU°ERHEATER    TUBE    METAL    TEMPERATURE  CS^C^190 

*  C  SMC 4  2 CO 
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MLE:     CS^P               FORTRAN       °1  NAVAL    ?3  ST^RAD  LA  TE    SCHOOL 

T*i.=  INTGRL<  T*k»C,GTWfc)  CSMC4210 

*  „_  CSvCf220 

*  CCfF'JTE    Th£     IMPLICIT     EQUATION     FOR    TlTil      STFAM  CS*C4230 

*  FLCh    RATE    f-RC  f     EDILtR  L    5,tA^                                             CSMC4240 

*  .  CS^04250 
M»*  =  l*FHNVMJ,vMNt,*MMi*p)  cSMC42to 

»  CSV042  70 

*  COMPUTE    THE    SJPFRhEATEP    OUTLET  P*ESSU3E                                                             CSMC4280 

*  CSMC4290 
£>MN=SO»T  (  OS4T  •*2.0-  2.0*KTN2*K0N1**<'MM4*2«  0)  CSM0<t300 

nwm  1 1  =  pm**kcn4*  valve  c  sw  C43io 

mmmi  l=Mf\jr.jo  CSM04320 

ppk]  *.F=MV«II+MVMHI  CSM04330 

*  CSM04340 

*  EQUATE    THE    FLUE    GAS    MASS    FLOW    RATES  .    CSM04350 

*  CSNC4360 
"RR=MOQ  CSM04370 
mSs=mrr  cS  ^04360 
NTT=*SS  CSMC4290 
f^LL=MTT  CSMG44GQ 

*  C5MQ4410 

*  SET    hiScR    INLET    DENSITY    EQUAL    TO  CSMC442G 

*  SATLRATED    LIQJIO    DENSITY  CSM04430 

*  CSVC4440 
FHCEE  =  RHCF  CSMC4450 
9H0JJ=RHCF  CS"04460 
CHruC=  RhOEE  CSMU4470 
RHCKh  =  RHCJJ  CSMC4480 

*  •  CSM04490 

*  ECU^te    CCWNCC^bR    FLr«     RATES    TO  °ISE°    EMT^A'^.E    FLOW                              CSMC4500 

*  RATES    AND    SET     RISER    ENTRANCE    FLCW    ?AT£S    TJ    THE  CSMG4510 

*  FLOW    RATES    AT    THE     INITIAL     THE  CS*04520 

*  .  CSMCA530 
VHH=MGG  CSM04  540 
*CC=VDC  CSWC4550 
►GG=MJ  CS'-'C4560 
M3D=MFFC  CS"i04570 
*JJ=wLLO  CSM04580 

*  •  CSUC45S0 

*  C3M=>UTE    THE    DERIVATIVE     JF    AVERAGE    RISER    DENSITY  CSMQ46G0 

*  C  S  M  C  4  1 1 0 

*  CSMC4620 

*  COMPUTE    THE    STEAM    MASS    RATE    EQUATION    FOR    THt    STEAM  CS^C4630 

*  DR'JN  CS,''C4fc4C 
DDSTM=*FF*VFF  +  XLL*,'LL-,,,CCND-MMM  CS  v 04650 

*  CS^- 04660 

*  CC^FUTE    THE    CRU*    STEA*    MASS  CSMC4670 

*  CS  '^04680 
CSTM  =  INTCR I (C^T*0,D0STV  J  CSVC46S0 

*  CSM  0*7  GO 

*  C3MPIJTE    THE    VOLUME    3F    STEAM    IN  THE     STEAM    ORjx                                             CS"G4710 

*  C SMC 4 7 20 
vn|_STV=VQLCRM-QMASL/RHOF                                        ..  CSV04730 

*  CSJO4740 

*  CCMPLTE    THE    DENSITY     OF    STEAM    IN    THE    STEAM     C^UM  C5MC4750 

*  CSN04760 
RHCST^  =  CcO/VCLSTM  CSMC4770 
RHOMM=RHCSTf  CS'-'C47£0 
W«M  =  1  .Q/RHCMM  CS"04790 
VFGM*=VV^-VF  CS^C48C0 

*  CSMC4810 

*  CCMOUTE    THE    STEAM    Dc UM    STEAM    OUTLET    PRESSURE  CS"C4820 

*  CSM04830 
"°ESSW=5  24.C/( VFGMM+. 1)  CS "048  40 
PRCCEDURE  F^N  =  FILTF5  (FKESSM)  CSWC4850 
lFlPt>FSSM.LT.G.Q}CALL  CE3UG  (3,0.0)  CSM04860 
FMw=FFESSM  CSVC4870 
Ef\JPPCCEClPE  CSMC4880 

*  CSM04890 

*  COMPUTE    THE    STEAM    ORJM    STEA"    3UTLET    TEMPERATURE  CSMG4900 
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FILE:    CS^IP 


FQRTRAM       PI 


NAVAL    PUSTGRAOlATE    SCHCCl 


TM*=cXD(  . 22  1 5 1*AL3G(P'4M)<.4.  771 23) 


SOLVE    FOR    THE    ECCNQMI2ER    FEED    CUTLET    ENTHALPY 


FEECRATE    EQUATION 


TWO    PHASE    PLC*    MULTIPLIER    EQUATIONS 


* 

* 

HBE=(KEB*HAA3+C8)/MBE 
* 
4 
* 
*BB=M830 

* 

PiiPAVE=2  4.7Q**XFF**2  .J  -6  .  5<Jtd*XF  F+  .  9  77b 
RJrA  VK=2A.  79^*  XL L**  2.0- 6.  50  6  6* >L L+ . 9  776 
c  ACL  t=L3.45fct*xFF**^.vJ«-io.  4  9  44*XFF-.  C0Q0  7 
RACL  J  =  15.<+56-+*XLL**2.J  +  13  .494*  *X  LL-.00  00  7 
P.FRICE*-  34.  C £22*  XFF** 2.  C+ 2 3.  7164*XFF +.8734 
PFPiCK  =  -.i<*.0a22*XLL**<i.C*2j.716<t*XLL+.a734 

*  STATE     POINT    EOLATIONS 

* 

nSAT=OH 

Hr=riSAT 

bF  L=  9  22.  15-0. 435 1 6*PSA  T  *  1.  717  E-04*°S  AT  **2.C- 4.  2 19E- 0  6... 

*PSAT**3.0 
|-V  =  HS*T*HFG 
HH  CF  =63.  8 -0.31 73  1*TS  AT  4l  .1  32  E-0  5  *TS  AT**2  .0-6  .  786E-03  .  .  . 


CSMC4910 

CSMCK920 

CS*C4930 

CSM  C4  9  40 

CS  "C4950 

CSMC4960- 

CSMC4970 

csvg49bo 

CS^C4990 
CS  ""0  50  00 
CS^CbOlO 
CS'*C5C20 
CSM3  50  3O 
CS"C6C40 
CS"Q5050 
CS  ^05060 
CSMC5070 
CS  "105080 
CSV0509O 
C  SMC 5  100 
CS VC5110 
CSWQ5120 
CSMC5130 
CS"05140 
CSVC5I50 
CS"'05l60 
CS  *C5170 
CSMC5180 
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APPENDIX   C 


i  .    Anticipated  Performance 


Rai  of  Operation  -Per  Cent 
I  oil  Steam  Generated  lb/hr 
jjrheated  Steam  Ibs/hr. 
L;a:perlieated  Steam  Ibs/hr. 
Boi  r  Drum  Pressure  psig 
jujrhcater  Outlet  Pressure  psig 
jujrheater  Outlet  Temperature  °F 
[  e:perheated  Steam  Outlet  Pressure  psig 
Deipcrheated  Steam  Outlet  Temperature  °F 
Eccomizer  Inlet  Pressure  psig 
Eccomizer  Inlet  Temperature  °F 
Eccomizer  Outlet  Pressure  psig 
Liccomizer  Outlet  Temperatuie  °F 
Jjsig  Air  Inlet  Temperature  °F 
Toil  Air  Flow  lb/hr 
foil  Oil  Flow  lb/hr 
\n  ;ipated  Efficiency  % 
^aanteed  Efficiency  % 
Kantion  snd  Unaccounted  for  Losses  % 
tx<ss  Air% 
Cann  Dioxide  "c 
NJiber  of  Burners  in  Operation 
ihittle  or  Non-Throttle  of  Air  Doors 

Dra  Loss  -Total  Inches  Water 
rnijgh  Double  Casing 
Hirjgh  Burner  Register 
fiirugh  Boiler  &  Superheater 
Flu  jgh  Economizer 

jas'emperature  Leaving  Superheater  Screen  °F 
jas'emperature  Leaving  Superheater  °F 
jas'emperature  Leaving  Main  Bank  °F     ' 
jas'emperature  Leaving  Economizer  °F 
lea Release  KB/Hr/Sq.  Ft.  Radiant  Heat  Absorbing 
Surce 

lea Release  KB/Hr/Sq.  Ft.  Total  Heating  Surface 
-a  Release  KB/Hr/Cu.  Ft.  Furnace  Volume 
:urice  Heat  Absorption  KB/IIr./Sq.  Ft 
iea  Absorption  First  Water  Scieen  Row  KB/Hr/Sq. 

Ft.  (Max.) 
!'a Absorption  Maximum  KB/IIr./Sq.  Ft 

Dr;  Losses  (Full  Power)  Based  On  241  Cu.  Ft.  of  100°F  Air/Lb.  of  Fuel  Oil 
*Dift  Losses  (Overload)  Based  On  24!  Cu.  Ft.  of  I0C°F  Air/Lb.  of  Fuel  Oil 
**l  aft  Losses  (Overload)  Based  on  260  Cu.  Ft.  of  68°F  Air/Lb.  of  Fuel  Oil 


Design 

Rated 

Maximum 

Endurance 

Full 

Intermittent 

Endurance          Boiler 

2  Boilers 

Power 

Power 

1  Boiler 

Overload 

56 

83 

97 

114 

120 

195,000 

290,000 

340,000 

400.000 

420.000 

155,000 

235,000 

265.000 

300,000 

320.000 

40,000 

55,000 

75,000 

100,000 

100.000 

690 

690 

690 

690 

690 

665 

030 

610 

580 

570 

912 

914 

980 

898 

895 

651 

606 

558 

480 

470 

635 

659 

680 

690 

690 

715 

726 

733 

749 

754 

250 

250 

250 

250 

250 

703 

705 

703 

708 

710 

355 

372 

372 

377 

380 

100 

100 

100 

100 

100 

242,229 

368.092 

433,589 

513,154 

540,952 

14^33 

21,781 

25,656 

30,364 

32,009 

85.7 

84.4 

83.8 

83.0 

82.7 

85.7 

84.4 

83.8 

83.0 

82.7 

.99 

.99 

.95 

1.01 

.98 

15.0 

15.0 

15.0 

15.0 

15.0 

13.0 

13.0 

13.0 

13.0 

13.0 

6 

6 

6 

6 

6 

N.T. 

N.T. 

N.T. 

N.T. 

N.T. 

13.74 

33.98* 

47.96 

69.51 

76.22**93.48 

1.33 

3.09 

4.28 

6.00 

6.67        10.22 

2.60 

6.00 

8.50 

11.50 

13.00      15.88 

4.83 

12.29 

18.14 

25.61 

28.15     38.32 

4.98 

12.60 

17.04 

26.40 

28.40     29.06 

2471 

2594 

2642 

2686 

2699 

1826 

1961 

202  J 

2085 

2103 

694 

774 

812 

856 

872 

373 

427 

453 

484 

498 

447 

680 

800.6 

948 

999 

15.5 

23.6 

28.0 

32.9 

34.7 

183.9 

279.5 

329.2 

389.6 

410.8 

128.4 

171.0 

191.1 

213.19 

221.5 

210 

231 

255.3 

288 

256 

210 


231 


255.3 


288 


256 
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II.      Anticipated  Metal  Temperature  Degree  F. 


Water  Screen  Tubes  Outside 
Water  Screen  Tubes  Inside 
Superheater  Tubes  Outside  -  Maximum 
Superheater  Tubes  Inside  -  Maximum 
Maximum  Inner  Casing  Temperature 
Maximum  Outer  Casing  Temperature 

at  way  of  Structural  ties 
Maximum  Outer  Casing  Temperature 


Design 

Maximum 

Endurance 

Rated 

Intermittent 

Endurance 

Boilei 

2  Boilers 

Full  Power 

Power 

1  Boiler 

Overload 

644 

656 

661 

667 

680 

554 

554 

554 

554 

555 

1022 

1046 

1049 

1048 

1048 

1044 

1020 

1019 

1014 

1012 

SOO 

825 

835 

850 

870 

350 

350 

350 

350 

350 

145 

145 

145 

145 

145 

Tube  Data 


O.D. 

M.W.T. 

No. 

Side  Wall  and  Roof 

2" 

.134" 

71 

Rear  Wall 

2" 

.134" 

54 

Front  Wall 

2" 

.134" 

22 

Screen  Bank 

2" 

.134" 

102 

Superheater 

1.5" 

.120 

22 

Main  Bank 

2" 

.134" 

34 

Main  Bank 

1" 

.085" 

2808 

Economizer 

2" 

.180" 

182 

Desuperheater  (in  Water  Drum) 

2" 

.220" 

6 

Downcomers 

8  5/8" 

.483" 

1 

» 

Downcomers 

10  3/4" 

.519" 

5 

Downcomers 

12  3/4" 

.601 

2 

Risers 


.500" 


Material 

M1L-T-16286CL.  A* 
MIL-T-16286CL;  A* 
MIL-T-16286CL.  A* 
MIL-T-16286CL.  A* 
MIL-T-16286CL.  E 

MIL-T-16286CL.  ,v* 
MIL-T-16286CL.  A* 
MIL-T-16286CL.  a 
SA-268  TP430 


Schedule  80  Pipe 
ASMESA-106-B 

Schedule  80  Pipe 
ASMESA-106-B 

Schedule  80  Pipe 
ASMESA-lOd-B 

ASMESA-106-B 


These  tubes  may  be  MIL-T-16286,  Class  A  (Seamless)  oi  MIL-T-1 7188  (Seamed)  electric  resistance  welded. 
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